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Making a combined use of bosonization and fermionization techniques, we build nonlocal trans¬ 
formations between dual fermion operators, describing junctions of strongly interacting spinful one¬ 
dimensional quantum wires. Our approach allows for trading strongly interacting (in the original 
coordinates) fermionic Hamiltonians for weakly interacting (in the dual coordinates) ones. It enables 
us to generalize to the strongly interacting regime the fermionic renormalization group approach to 
weakly interacting junctions. As a result, on one hand, we are able to pertinently complement the 
information about the phase diagram of the junction obtained within bosonization approach; on the 
other hand, we map out the full crossover of the conductance tensors between any two fixed points 
in the phase diagram connected by a renormalization group trajectory. 

PACS numbers: 71.10.Pm, 72.10.-d, 73.63.Nm 


I. INTRODUCTION 

Many-electron systems, such as conduction electrons in a metal, are typically well-described by Landau’s Fermi 
liquid theory (LFLT), which is capable of successfully accounting for the main effects of electronic interaction-lr— . 
LFLT’s basic assumption is that, close to the Fermi surface, low-energy elementary ’’quasiparticle” excitations of the 
interacting electron liquid are in one-to-one correspondence with particle- and hole-excitations in the noninteracting 
Fermi gas. This corresponds to a nonzero overlap between the quasiparticle wave function in the Fermi liquid and the 
electron/hole wave function in the noninteracting Fermi gas, which is reflected by a coefficient Zp of the quasiparticle 
peak at the Fermi surface in the spectral density of states that is finite though, in general, < 1 (Zp = 1 corresponds 
to the noninteracting limit). The ’’adiabatic deformability” of quasiparticles to electrons and/or holes by smoothly 
switching off the interaction allows, for instance, to address transport in a Fermi liquid in a similar way to what is 
done using scattering approach in a noninteracting Fermi gas, etc£. 

LFLT is grounded on the possibility of obtaining reliable results by perturbatively treating the electronic 
interaction^. This is strictly related to the small rate of multi-particle inelastic processes, in which, due to the 
interaction, an electron/hole emits electron-hole pairs. While this is typically the case in systems with spatial dimen¬ 
sion d higher than one, in one-dimensional systems, the proliferation of particle-hole pair emission at low energies 
leads to a diverging corresponding rate, which makes the quasiparticle peak disappear [Zp = 0). As a result, the 
interaction cannot be dealt with perturbatively and one has rather to resort to nonperturbative techniques, allowing 
for summing over infinite sets of diagrams^ - —. The breakdown of LFLT in one-dimensional interacting electronic 
systems (interacting ’’quantum wires” - QW’s) reveals itself in a series of physical contexts, the more remarkable of 
which is the transport across an interacting QW in the presence of a constriction, or a weak link, or the conductance 
of a junction of QWs, with or without spinier- 3 . Moreover, due to the remarkable correspondence between fermionic- 
and spin-bosonic one-dimensional lattice models encoded in the Jordan-Wigner fermionic representation of quantum 
spin-1/2 spin operators^, a number of bosonic realizations of one-dimensional models of interacting fermions have 
been proposed in bosonic systems, such as quantum spin chains^r— , quantum Josephson junction networks^ - —, or 
topological Kondo-type systems^— . 

Formally, interacting electrons in one dimension are commonly treated within Tomonaga-Luttinger liquid (TLL)- 
approac h 39 : 40 . TLL formalism provides a general description of low-energy physics of a one-dimensional interacting 
electronic system in terms of collective bosonic excitations (charge- and/or spin-plasmons). Electronic operators are 
realized as nonlinear vertex operators of the bosonic field s™^ . As for what concerns transport properties, the most 
striking prediction of the TLL-approach is possibly the power-law dependence of the conductance on the low-energy 
reference scale (’’infrared cutoff”), which is typically identified with the (Boltzmann constant times the) temperature, 
or with the (Fermi velocity times the) inverse system length (’’finite-size gap”) in a dc transport measurement, or 
with eV, with V being an applied voltage, in a nonequilibrium experiment. While TLL-formalism poses no particular 
constraints on the strength of the ’’bulk” interaction within the quantum wires, it suffers of limitations, when used 
to describe transport across impurities in an interacting quantum wire, or conduction properties at a junction of 
quantum wires (generically referred to, in the following of this section, as quantum impurities in the quantum wire). 
Parametrizing the spatial coordinate in each wire with a:(> 0), within the TLL-framework, a junction of quantum wires 
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is mapped onto a model of AT-one dimensional TLLs (one per each wire), interacting with each other by means of a 
’’boundary interaction” localized at x = 0. Dealing with such a class of boundary problems requires pertinently setting 
the boundary conditions on the plasmon fields at x = 0. While in some very special cases the boundary conditions 
can be written as simple linear relations between the plasmon fields, in general they cannot. This is a consequence 
of the nonlinearity of the relations between bosonic and fermionic fields: even linear conditions among the fermionic 
fields are traded for highly nonlinear conditions in the bosonic fields. As a consequence, except at the fixed points of 
the boundary phase diagram, where the boundary conditions are ’’conformal”, that is, linear in the bosonic fields, the 
boundary interaction can only be dealt with perturbatively, with respect to the closest conformal fixed poin t 43 : 44 . With 
very few remarkable exceptions, where the junction model Hamiltonian maps onto some exactly solvable model^c— , 
the only way for recovering ’’global” (i.e., not necessarily in the vicinity of a conformal fixed point) informations about 
the phase diagram and the corresponding scaling of the conductance, is by just making educated guesses from the 
global topology of the fixed-point manifold of the phase diagram. Recently, the bosonization approach in combination 
with zero-temperature numerics has successfully been employed to compute the junction conductance by relating it 
to the asymptotic behavior of certain static correlation functions. Correspondingly, the length scale over which the 
asymptotic behavior emerges has been worked out, as wel l 48 : 49 . Also, numerical calculations of the finite-temperature 
junction conductance can be performed by using quantum Monte Carlo approach^ or, likely, by implementing some 
pertinently adapted version of the finite-temperature density matrix renormalization group approach to quantum spin 
chain s 51 : 52 . 

Alternatively, one may not be required to give up using fermionic coordinates, by employing a systematic renormal¬ 
ization group (RG) procedure to treat the effects of the bulk interaction on the scattering amplitudes at the junction. 
Among the various possible ways of implementing RG for junctions of interacting QWs, the two most effective (and 
widely used) ones are certainly the poor man’s fermionic renormalization group (FRG)-approach, based upon a sys¬ 
tematic summation of the leading-log divergences of the 5-matrix elements at the Fermi momentum and typically 
yielding equations that can be analytically treated ^ 15 : 18 , and the functional renormalization group (fRG)-approach, 
based on the functional renormalization group method, leading to a set of coupled differential equations for the vertex 
part that, typically, can only be numerically treated^ - — . Both approaches are expected to apply only for a sufficiently 
weak electronic interaction in the quantum wires and, in this sense, they are less general than the TLL-approach, 
which applies even for a strong bulk interaction. Nevertheless, at variance with the TLL-approach, a RG-approach 
based on the use of fermionic coordinates leads to equations for the 5-matrix elements valid at any scale and, thus, 
it allows for recovering the full scaling of the conductance, all the way down to the infrared cutoff. 

Aside from the remarkable merit of providing analytically tractable RG-flow equations, the FRG-approach, when ap¬ 
plied to interacting spinful electrons, also accounts for the backscattering bulk interaction, which is usually neglected 
in the TLL-framewor k 13 : 14 . Moreover, it can be readily generalized to describe junctions involving superconduct¬ 
ing contacts, at the price of doubling the set of degrees of freedom, to treat particle- and hole-excitations on the 
same footing^. Nevertheless, FRG is known to suffer of limitations, when used to describe the crossover towards 
’’healed” fixed points, where quantum impurities are renormalized towards boundary conditions corresponding to per¬ 
fect conduction properties. In particular, this leads to a scaling equation for the conductance which, in the pertinent 
asymptotic regime, is not consistent with the one obtained from the bosonization approach. Using the fRG-approach 
allows for taking care of this flaw, as the fRG-technique typically takes into account the mutual feedback from all the 
single-fermion scattering channels, and not only those from scattering processes between different Fermi points (see 
Ref. [HI for a detailed discussion of this point and for a careful comparison between FRG- and fRG-approaches). In 
fact, fRG appears to be generically more accurate than FRG (which can be in fact recovered from fRG under suitable 
approximations^). Both approaches suffer, however, of the limitation on the bulk interaction, which must be weak, 
in order for the technique to be reliably applicable. 

To overcome such a limitation, in this paper, we study a junction of spinful interacting QWs by making a combined 
use of bosonization and fermionization, that is, we go back and forth from fermionic to bosonic coordinates, and vice 
versa, to build ” dual-fermion” representations of the junction in strongly interacting regimes. In resorting from a 
bosonic to a fermionic problem, our approach is reminiscent of the refermionization scheme used in Ref. [64j to discuss 
the large-distance behavior of the classical sine-Gordon model at the commensurate-incommensurate phase transition. 
Specifically, in^ the refermionization allows for singling out at criticality the low-energy two-fermion excitations from 
the one-fermion ones and to prove that the latter ones keep gapped along the phase transitions and do not contribute 
to the large-distance scaling of the correlations. At variance, in our case it is the second of a two-step process, that 
ends up again into a fermionic ” dual fermionic” model for the strongly-interacting system. The guideline to construct 
the appropriate novel fermionic degrees of freedom is to eventually rewrite the relevant boundary interactions as 
bilinear functionals of the fermionic fields. Specifically, moving from the original fermionic coordinates to the TLL- 
bosonic description of the junction, we are able to warp from the weakly interacting regime to different strongly 
interacting regimes. Therefore, at appropriate values of the interaction-dependent Luttinger parameters, we move 
back from the bosonic- to pertinent dual-fermionic coordinates, chosen so that the relevant boundary interactions 



3 


are bilinear functionals of the fermionic fields. Our mapping between dual coordinates is actually preliminary to the 
implementation of the RG-approach. Therefore, in principle, it could be equally well applied to extend both FRG 
and fRG to the strongly-interacting regime. Since, however, fRG usually requires resorting back to lattice models to 
implement a systematic numerical treatment of the flow equations for the interaction vertices, which goes beyond the 
scope of this workaa—, we rather prefer to complement our dual mappings with a pertinently adapted version of the 
FRG-approach. 

The RG-approach formulated in fermionic coordinates, such as FRG, suffers of the limitation that it requires that 
relevant scattering processes at the junction are fully encoded in terms of a single-particle 5-matrix. While this is 
certainly the case at weak bulk interaction, a strong attractive interaction in either charge-, or spin-channel (or in 
both) is known to stabilize phases (RG attractive fixed points) at which two-particle scattering is the most relevant 
process at the junction— ~ 12 ’ 20 ’ 21 . Just because of the way it is formulated, the FRG-approach fails to describe many- 
particle scattering processes, even after improvements of the technique that allow to circumvent the constraint of 
having a small bulk interaction^—. Resorting to the appropriate dual-fermion basis allows us to describe within 
the FRG-approach also fixed points stabilized by many-particle scattering processes, as well as fixed points whose 
properties have not been mapped out within the TLL-framework in terms of a rotation matrix such as, for instance, 
the mysterious-fixed point in the three-wire junction of spinless quantum wires studied in Ref. [2l| and its counterpart 
in the junction of spinful quantum wires. Moreover, in computing the conductance tensor along the RG-trajectories 
connecting fixed points of the phase diagram, we show how our approach, while being consistent with the TLL- 
approach in the range of parameters where both of them apply, on the other hand allows for complementing the 
results of Refs. GMIH about the two-wire and the three-wire junction, with a number of additional results about 
the topology of their phase diagram and their conductance properties. 

The paper is organized as follows: 

• In Sec. uu we apply our duality-complemented FRG-approach to a junction of two interacting spinful quantum 
wire, discussing the results and comparing them to those obtained within the bosonization approach^—; 

• In Sec. m we apply our approach to a junction of three interacting spinful quantum wire. We first compare 
our results with those obtained within the bosonization approach^ and, therefore, we show how our technique 
allows for mapping out the full crossover of the conductance at the junction even in strongly-interacting regions, 
typically not accessible with a fermionic approach; 

• We summarize our results and discuss possible further developments of our research in Sec. EY1 dedicated to 
the concluding remarks of our paper; 

• In the various appendixes we review mathematical techniques that are crucial for our derivation. Specifically, 
in Appendix [A] we review the derivation of the FRG-equations for the 5-matrix, in Appendix [B] we review the 
basic bosonization rules for interacting one-dimensional quantum wires; in Appendix [U] we provide some basic 
elements of linear transport theory for junctions of one-dimensional quantum wires. 


II. DUAL FERMIONIC VARIABLES AND RENORMALIZATION GROUP APPROACH TO THE 
CALCULATION OF THE CONDUCTANCE AT A JUNCTION OF TWO SPINFUL INTERACTING 

QUANTUM WIRES 


To introduce and check the validity of our approach, in this section we discuss a junction of two interacting spinful 
quantum wires. This appears to be quite an appropriate place to test our technique: indeed, the two-wire junction has 
widely been studied in the past, both within the bosonization approach^—, and by means of standard RG techniques 
for a weak bulk interaction, either using the FRG-approach ^ 15 ’ 18 , or the fRG-method£2r— . The two-wire junction 
of spinful quantum wires is described by the (K = 2) bulk Hamiltonian iLeuik — Ho + H mt , with 


Hq = —iv 


dx R,jA x ) 

j =i Jo 




L,j,cr 


C x)d x ^Lj, a (x ) j 


(1) 


with L being the wire length (eventually sent to infinity at the end of the calculations) and the interaction Hamiltonian 
given by 
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The various interaction strengths appearing in Eq. 0 are defined as 


ST i,|| = v j,n( 2k p) = v j,u( 2k p) 

9j, 2,|| = ^j,tt(°) = ^,44 (°) 
flj.i.-L = v j,n( 2k p) = v jm ( 2k p) 

9i, 2,± = ^,n(0) = Vat( 0 ) > ( 3 ) 

with Vcr tC r'(k) being the Fourier modes of the two-body interaction ’’bulk” interaction potential in the quantum wires 
(see Appendix [XI for the derivation and discussion of Eqs. (I2I3D .1 For a weak bulk interaction, the most relevant 
contribution to Hb is given by a linear combination of the operators defined as 


B (j,j'),cr,(x,x'){ 0) = VE^WV’X'J',^ 0 ) 


( 4 ) 


with X , X' = L 1 R. Assuming equivalence between the two wires and a spin-symmetric and spin-conserving boundary 
interaction, Hb can be generically written as 


H b 


E 

X,X'=L,R 


E {[ r A',X' J B(l,2),<T,(X,.Y')(0) + h.C.] + ^2 VX,X' B <J,j),<T,(X,X'){ 0)} 

cr 3 = 1,2 


( 5 ) 


In addition to the contributions reported in Eq. 0, terms that are quadratic (or of higher order) in the B’s can in 
principle arise along RG-procedure, even if they are not present in the ’’bare” Hamiltonian. For instance, the simplest 
higher-order boundary interaction terms consistent with spin conservation at the junction, H 2 0 and Hq 2 , are given 

byi2-i£ 


#2,0 — ^2,0 ^2 {-®(X,X'),t,(l,2)(0)-S(X',X)4,(l,2)(0) + h.c.} 

X,X'=R,L 

Ho ,2 = Vo ,2 ^2 {-®(A,X'),t,(l, 2 )( 0 )S(X',X) 4 ,( 2 ,l)( 0 ) + h.c.} . (6) 

X,X'—R,L 

As it can be shown using the bosonization approach, for a weak bulk interaction, higher-order operators such as 
those in Eqs. 0 are highly irrelevant operators and, accordingly they are typically ignored and one uses for Hb the 
formula in Eq. 0- Physically, this means that the relevant scattering processes at the junction consist only of one 
single particle/hole scattered into one single particle/hole, such as those drawn in Fig. |T| (a). These processes are 
fully described by the single-particle S-matrix, for which the renormalization group equations can be fully recovered 
using the technique we review in Appendix [A] The symmetry requirements listed above imply that the single-particle 
S'-matrix takes the block-diagonal form 


'-’(3>);0',<r , )(&) — &o,o' , 


( 7 ) 


with the S(k )-matrix being given by 


S(k) 


Tk tk 
tk Tk 


( 8 ) 


and rfc and tk respectively corresponding to the amplitude for the particle to be backscattered in the same wire, 
or transmitted into the other wire. In the following we will pose no particular constraints on the ry.’s and the f/.’s, 
except that, near the Fermi points, they are quite flat functions of fc, without displaying particular features, such as 
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a resonant behavior: accordingly, we assume that the amplitudes are all computed at the Fermi level and drop the k, 
label (this is a specific case of the general assumptions on the behavior of the S'-matrix elements near by the Fermi 
surface that we make in Appendix [Aj . To write the RG-equations for the S'-matrix elements, one needs the Friedel 
matrix F— which, in this specific case, is given by 



" f3r 
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(9) 


with P ■= —( 9i\\ ~ 9i± + < 72 11) • Taking into account the symmetries of the S- and of the F-matrix, the RG-equations 
for the amplitudes r, t are obtained in the form 

^ = f ( r - r M 2 - r *t 2 ) = P r \tf 

= -Pt\r\ 2 = -p(t-t\t\ 2 ) . (10) 

Equations m must be supplemented with the RG-equation for the running strength /3, which is given by 

= + 29l ' ± ^ 92 ' ± ~ 92 ’ ll+91 ^ ■ (U) 

(See Eqs. (12131) for the definition of the bulk interaction strengths <7y,i(2),_L-> <7y,i(2),|| : here we drop the wire index j as 
the interaction strengths are assumed to be the same in each wire.) Equations COD, together with Eq. m and Eqs. 
(1A10I) for the running interaction strengths, constitute a closed set of equations, whose solution yields the scaling 
functions r(D),t(D). From the explicit formulas for the running scattering amplitudes, one may readily compute 
the charge- and the spin-conductance tensors, using the formulas derived in Appendix 1C II As a result, due to the 
symmetries of the S’-matrix, the charge- and the spin-conductance tensor are equal to each other and both given by 


G C (D)=G S (D) 


£ ~ G(D) G{D) 
G{D) £ - G(D) 


( 12 ) 


with G{D) = ^\t(D)\ 2 . An explicit analytical formula can be provided for G(D) in some simple cases such as, for 
instance, if gi± is fine-tuned to 0. In this case, as it arises from Eq. CD. P keeps constant along the RG-trajectories 
and, therefore, one may exactly integrate Eqs. m for r(D) and t(D), obtaining 


r{D) 

t{D) 


t o \D/D o \- 0 

yJ\ro\ 2 + \t o \ 2 \D/D o \~ 20 

to |-P/A/ 

yf\ro \ 2 + \to\ 2 \D/D 0 \ 2 P 


(13) 


with ro,to corresponding to the ’’bare” scattering amplitudes in Eq. (j8]). Another case in which an explicit analytical 
solution can be provided corresponds to having gi_L = g±\\ = g\ and ^ 2 _l = < 72 1| — 92 - In this case, the set of Eqs. (1A10I) 
collapse onto a set of two equations for gi(D), g 2 {D) which can be readily integrated, yielding the running interaction 
strengths 


9i(D) = 


9 1 


1 + ^ In 
irv D 
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7 TV D 

and P(D) = [ 92 (D) — 2gi(D)\/(2nv). Once P(D) is known, Eqs. (fTOll can be integrated, yielding 
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FIG. 1: Sketch of possible scattering processes at a two-wire junction (note that incoming particles from wire j can either be 
scattered into the same wire, or into a different wire): 

(a) Single-particle/single-hole scattering processes. These are determined by Htun in Eq. © and are fully described in terms 
of a single-particle S'-matrix; 

(b) Many-body scattering processes in which one particle/one hole is scattered into two particles and one hole/two holes 
and one particle. These processes can be induced by boundary interaction Hamiltonians such as those in Eq. © and their 
proliferation requires resorting to a bosonic Luttinger-liquid description of the junction; 

(c) Scattering processes for a particle-particle and for a particle-hole pair. These are again determined by the Hamiltonians 
in Eq. © and are the only allowed processes in the presence of a strong repulsive (attractive) interaction in the spin (charge) 
channel, and vice versa. On pertinently defining new fermionic coordinates, they can still be described in terms of a "single-pair” 
S'-matrix. 


with 7 = (— 3 - + 32 ) /(27n>), Tq = |fo| 2 ,-Ro = |fo| 2 - As an example of typical scaling plots for G{D) in the simple 
cases discussed before, in Fig. Owe plot G(D)7 r/e 2 vs. \n(D 0 /D), as from Eq. (fTdl) (panel (a)) and from Eq. (fTSl) 
(panel (b)), with the values of the parameters reported in the caption. Consistently with the results obtained within 
Luttinger liquid frameworbilT— , G(D)7 r/e 2 either flows to 0 for an effectively repulsive interaction (/3,y > 0), or to 2 
(the maximum value allowed by unitarity), for an effectively attractive interaction (/?, 7 < 0). For general values of the 
interaction strengths, the equations have to be numerically integrated. In Fig. © we provide some examples of scaling 
of G(D) vs. \n(Do/D) in the general case. It is important to stress^ that, due to the nontrivial renormalization group 
flow of the interaction strengths, the flow of G(D) can be a nonmonotonic function of D for some specific values of the 
interaction strengths. It would be interesting to check such a feature in a real life two-wire junction: remarkably, this 
prediction is only obtained within the FRG-approach, in which it is possible to account for the flow of the running 
interaction strengths, as well. 

The possibility of mapping out the full crossover of the conductance as a function of the scale D is possibly the most 
important feature of the FRG approach. Yet, since, as we discuss to some extent in Appendix [Aj the validity of the 
FRG technique is grounded on the assumption that all the relevant scattering processes at the junction are described 
by the single-particle ff-matrfxpli T 15 ! 18 , it breaks down when attempting to recover the full crossover of the conductance 
towards fixed points where multi-particle scattering is the most relevant process at the junction, such as the strongly 
coupled fixed point stabilized by either # 2,0 or Hq 2 in Eq. Technically, what happens is that, as soon 

as boundary operators such as I? 2 ,o or Flo ,2 become relevant, the proliferation of low-energy many-body scattering 
processes such as the one we sketch in Fig. ffl(b) invalidates the single-particle S'-matrix description of the junction 
dynamics. Nevertheless, some many-body scattering processes can be strongly limited by having, for instance, a 
strong repulsive interaction among particles with the same spin and a strong attractive interaction among particles 
with the same charge. In the bosonic framework, this corresponds to having values of the Luttinger parameters in Eqs. 
(IB 1 lh such that g c >2,g s <C 1. Indeed, in this limit on one hand, the strong spin repulsive interaction forbids the 
single-particle processes described by Hb in Eq. (©) (at small values of the boundary coupling strengths tx,x' , Hx,x’ 
this can be readily seen from the explicit result for the scaling dimension of Hb computed within the bosonization 

approach, which is a 'B.weak = 1 — ^- 5 ^“, which becomes 1, corresponding to a largely irrelevant operator). On 

the other hand, one expects that the strong charge attraction stabilizes tunneling of composite objects carrying zero 
spin, such as two-particle pairs, as the one depicted at the left-hand panel of Fig. [I] (c). As a consequence, due to 
the fact that these are again one-into-one scattering processes, one expects that it is possible to choose the effective 
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ln(D 0 /D) 


ln(D 0 /D) 


FIG. 2: (a) Plot of G(D) vs. ln(Do/D) as from Eq. (1131) for To = 1 — Ro = 0.4, 0 = 0.35 (purple curve - corresponding to 
an effectively repulsive interaction), and 0 = —0.35 (blue curve - corresponding to an effectively attractive interaction), with 
To = \t 0 \ 2 , Ro = |to| 2 ; 

(b) Plot of G(D) vs. ln(Do/D) as given in Eq. (fT5T) for To = 1 — Ro = 0.4, gi/(2nv) = 0.2, 7 = 0.36 (purple curve - 
corresponding to an effectively repulsive interaction), and 0 = —0.36 (blue curve - corresponding to an effectively attractive 
interaction). 



b) 



FIG. 3: (a) Plot of G(D) vs. \n(Do/D), obtained by numerically integrating Eqs. (I10I11IA10I) for gqj (Do)/(2nv) = 

gi,±(D 0 )/(2ivv) = g 2 ,±{D 0 )/(2nv) = 0.3, g 2 \\(D 0 ) / (2irv) = 0, Ro = T 0 = 0.5; 

(b) Same as in panel (a), but with g 2 ,|| (0)/(27 tu) = 0.3, as well. 


low-energy degrees of freedom of the system to resort to a single-particle S'-matrix in the new coordinates. In fact, 
this is the idea behind the dual fermion approach we are going to discuss next. When resorting to dual fermion 
coordinates, an important issue is related to whether the vacuum states at a fixed particle number— for the original 
and the dual-fermions are the same. In fact, while dual fermion formalism only captures composite excitations e.g. 
two-particle states in the g c ~ 2,g s -C 1-regime, at such values of the parameters, these states are the only ones that 
at low energy are effectively able to tunnel across the junction (that is, the only ones whose tunneling is described 
by a non-irrelevant operator). So, as long as one is only concerned about states relevant for low-energy tunneling 
across the junction (that is, states relevant for the calculation of the dc-conductance tensor of the junction), one can 
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effectively assume that the fixed-particle number vacuum states are the same in terms of the new (dual) and of the 
old fermions. The definition of the dual fermion operators strongly depends on the boundary conditions of the various 
fields at the junction. Accordingly, in the following we define different dual fermion operators in different regimes of 
values of the boundary interaction, and eventually show that, whenever two different sets of dual coordinates apply 
to the same region, they yield the same results, as they are expected to. 

Let us begin with the weak boundary interaction regime. Referring to the bosonization formulas of Appendix |Bj 
this corresponds to assuming Neumann (Dirichlet) boundary conditions for all the $ (@)-fields in Appendix [Bl and, 
in addition, to equating the Klein factors so that = riR. a ,j, V<7. j. In this limit, one may respectively rewrite 

H 0 2 and I? 2 ,o i n bosonic coordinates as 


#2,0 = W 2 ,0 COs[$i jC (0) - $2,c(0)] 

#0,2 = ^0,2 COS[$1, S (0) - $2, s (0)] , (16) 

with V 2 ,o oc 14 ,0, Wo ,2 °c Vo, 2 - The scaling dimensions of the operators in Eqs. ca are respectively given by * 2,0 = 
2 /g c ,xo,2 = 2 /g s . Thus, in the regime of a strongly attractive interaction in the charge (spin)-channel and strongly 
repulsive interaction in the spin (charge)-channel, # 2.0 (# 0 , 2 ) may become the most relevant boundary operator at 
weak boundary coupling. The strategy of our dual fermion approach consists in defining a novel set of fermionic fields, 
in terms of which the operators in Eqs. (fTfil) are realized as bilinears, similar to the R-operators in Eq. (j4]). To be 
specific, let us introduce the center-of-mass and the relative fields in the charge- and in the spin-sector, respectively 
given by 


$ c<Y)(z) = ^=[$l, c ( S )(aO + $2,c( S )0*0] 
©c( S )(z) [®l,c(s) (*w) T ©2,c(s) (*w)] , 


(17) 


and 


Vc{s) (x) = -j=[$i, c ( s )(x) - $ 2 ,c( s )(a;)] 

$c(s){ x ) = [®l,c(s) ( x ) ~ 02,c(s)( :£ )] 

Next, let us perform the canonical transformation to a new set of bosonic fields, defined as 


(18) 


?c(s)0*0 


■V2 0 0 0' 


■ $ c(s) (x) ' 

®c(s) (**0 


0 

0 

0 


®c(s) ( x ) 

<£c( s )(z) 


0 0 V2 0 


Vc{s)( x ) 

^c(s) (*^) 


1 

0 

0 

0 

1_ 


- Vc(s)(x) - 


It is worth stressing that the transformation in Eqs. (1191) relate to each other bosonic operators at a given position in 
real space. Since the correspondence rules between the bosonic and the (original or dual) fermionic fields, summarized 
in Appendix iBl are local in real space, as well, one concludes that, written in terms of dual fermionic coordinates, 
the boundary interaction Hamiltonian Hr is still local and that the dynamics far from the junction can be fully 
encoded within dual fermion scattering states. Now, assuming g s -C 1, g c = 2 + Sg c , with \5g c \ -C 1, we see that 
g a -C 1 makes Uo ,2 strongly irrelevant. This fully suppresses spin transport across the junction and, therefore, we 
may just focus onto charge transport, ruled by # 2 , 0 - In fact, it appears that single-spinful particle-tunneling processes 
are already suppressed against two-particle pair tunneling processes as soon as g s < 2/3. As conservation of spin 
symmetry implies g s = 1, in order to realize the condition above one may, for instance, think of two coupled spinless 
interacting one-dimensional electronic systems (which could possibly realized as semiconducting quantum wires in the 
presence of spin-orbit and Zeeman interactions), with a mismatch in the Fermi momenta that prevents the interaction 
from opening a gap in the fermion spectrum. The two channels can, therefore, be regarded as the two opposite spin 
polarization, although without any symmetry implying g s = 1. To rewrite this latter operator as a bilinear functional 
of fermionic operators, we define the spinless chiral fermionic fields XR,ji x )>XL,j( x ) as 


XR ,j{x) = r] R j e i[*c(x)~(-iyM*)+&cW-(-iyMx)] 
Xj, j(x) = gL j 


( 20 ) 
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with j = 1,2 and with r]Rj->riLj being real fermionic Klein factors. ’’Inverting” the bosonization procedure outlined 
in Appendix [B] into a pertinent re-fermionization to spinless fermions, we find that the bulk Hamiltonian for the 
fermions is given by 


Hr- 


bulk 



dx {XR,j(x)d x XR,j(x ) - XL,j( x ) d =oXL,j(x)} 
'■ Xrj{x)xrj{x) :: XL,j( x )XL,j(x ) : , 


( 21 ) 


with the velocity u oc v. The Xi/fly-fields are the appropriate degrees of freedom to describe pair scattering at the 
junction in terms of a single-particle fj-matrix. In order to prove that it is so, we note that .& 2 ,o can be regarded as 
the bosonic expression for the boundary weak coupling limit of a tunnel Hamiltonian for the spinless fermions, H c tun , 
given by 


H c ,tun = «2,o{A t i(0)X2(0)+xk0)x 1 (0)} , (22) 

with Xi(0) = XR,j( 0) + Xlj(0)- While the strong spin repulsion sets the spin conductance tensor to 0, the charge 
conductance can nevertheless be different from zero, due to zero-spin pair-tunneling across the junction. Once the 
RG-flow for the S'-matrix elements describing x-fermion scattering at the junction has been derived as we did before, 
using the formulas we report in Appendix 1C II and the expression of the charge current operator in wire j in terms of 
the dual fermionic fields: 


J cj (x) = euV2{: Xrj( x )xrj( x ) ■ - ■ x[,j{ x )XL,j(x) :} 
we obtain that the charge conductance tensor scales according to 


with 


G C (D) 


i - G(D) G(D) 
G(D) G(D) 


G(D) 


7r 


T 0 \D/D 0 \-^ 

R 0 + T 0 \D/D 0 \-^ 


and the bare reflection and transmission coefficients respectively given by 


Ro 

T 0 


u 2 - v lo 

u 2 + 


2UV2A 


u 2,0 


(23) 


(24) 


(25) 


(26) 


In Fig. 01 we plot G(D) versus lii(I?o/^) in two paradigmatic cases, respectively corresponding to Sg c > 0 and 
to 8 g c < 0. To our knowledge, this is the first example of a full scaling plot of the conductance for a junction of 
strongly interacting one-dimensional quantum wires. While, on one hand, this shows the effectiveness of our approach 
in describing the crossover of the conductance towards the spin-insulating charge-conducting fixed point, on the other 
hand, one has also to prove the consistency of an effective theory strongly relying on the weak boundary coupling 
assumption with an RG-flow taking the system all the way down to the perfectly charge-conducting fixed point, 
corresponding to the strongly interacting limit of the boundary interaction^— . When 8 g c > 0, the relevance of U 2.0 
drives the system towards the strongly boundary interaction limit in the charge channel, corresponding to pinning 
<p c (0) and, accordingly, to imposing Neumann boundary conditions on "dc(O). Since $ c (0) does not appear in the 
boundary interaction, one assumes that it still obeys Neumann boundary conditions and, accordingly, that O c (0) is 
pinned at a constant value. We now prove that these boundary conditions are recovered by taking the strongly coupled 
limit of Ff Cj tim hr Eq. (1221) and using the refermionization rules in Eqs. (l20l) . Indeed, on making the strong-coupling 
assumption, |-2p| ~S> 1, as from Eqs. m 1 one obtains Rq —> 0,7b —> 1, that is, the boundary conditions correspond 
to perfect transmission from wire-1 to wire-2, and vice versa. In terms of the dual fermionic fields, this corresponds 
to the conditions 
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G(D) 7c/e 



FIG. 4: Plot of G C {D ) iis. ln(Do/.D) as from Eq. (1251) for To = Ro = 0.5 and for 6 g c respectively equal to 0.6 (blue curve) and 
to -0.6 (purple curve). 


XrM x ) = e * A 

Xr,i( x ) = e ~ iX Xl,2(~x) , (27) 

with A being some nonuniversal phase. From Eqs. one sees that Eqs. (1271) imply Dirichlet boundary conditions 

at x = 0 for both (p c (x) and 0 c (x), with the dual fields $ c (a;), $ c (x) obeying Neumann boundary conditions. This is 
exactly the same result one would obtain working in bosonic variables by sending to oo the interaction strength V 2 ,o in 
Eq. (flbl) . Due to the strong repulsion in the spin channel, such a fixed point corresponds to perfect transmission in the 
charge channel, but perfect reflection in the spin channel, that is, it must be identified with the non-symmetric charge¬ 
conducting spin-insulating phase of Refs. [10H12I|. To conclude the consistency check, we note that, on alledging for 
additional backscattering contributions to H c j U n of the generic form /j-ix'r i(0)X/? ^0) + 1*2x 1, i(0)xl (0) (which play 
no role at weak coupling) and using again Eqs. (1271) . one obtains the bosonic operators 

#c,tun ~ ftcos [^c(0)] , ( 28 ) 

with /i being some nonuniversal constant. Equation (1281) corresponds to the bosonic version of the leading boundary 
perturbation at the non-symmetric charge-conducting spin-insulating fixed point^ - — . 

Our approach also allows for analyzing the complementary situation in which j c <C 1 and g s ~ 2. In this case, 
one expects that the strong repulsion in the charge channel and the strong attraction in the spin channel stabilize 
single-pair tunneling processes at the junction such as those sketched at the right-hand panel of Fig. |T| (c), that 
is, tunneling of particle-hole pairs, with total spin 1. Again, for g s = 2, the S'-matrix describes single-particle into 
single-particle scattering processes, once it is written in the appropriate basis. To select the pertinent degrees of 
freedom, we therefore repeat the refermionization procedure in Eq. by just exchanging the charge- and the 

spin-sector with each other. Of course, charge- and spin-conductance are exchanged with each other, compared to 
the previous situation and, accordingly, the flow will be towards the charge-insulating spin-conducting fixed point of 
Refs. jloUl^ . An important remark, however, concerns the effects of a possible residual interaction, which, as we did 
before, can be in principle introduced for accounting for g s slightly different from 2. Indeed, a term in the ’’residual” 
bulk interaction Hamiltonian such as the one oc gyi,j_ in Eq. ([2]) , once expressed in terms of the fermionic fields in 
Eqs. (TSUI) would take the form 



with Tnj oc < 7 y, l ,_l , which would open a bulk gap in the single-x fermion spectrum, thus making the whole system behave 
as a bulk spin insulator. Therefore, in order to recover the correct physics of the charge-insulating spin-conducting 
fixed point, we must assume that all the <7j,i,_L are tuned to zero, which is typically the case when resorting to the 
bosonic approach to spinful electron&ilF— . 
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FIG. 5: Sketch of a three-wire junction of spinful quantum wires pierced by a magnetic flux <F. 


As we have just shown, resorting to pertinent dual-fermion operators allows for mapping out the full crossover with 
the appropriate energy scale of the charge and/or spin conductance of a junction in region of values of the interaction 
parameters in which one is typically forbidden to use the standard weak-coupling formulation of either FRG, or fRG. 
In the following, we apply our technique to the spinful three-wire junction studied using the bosonization approach 
in Ref. [23|, and will recover the full crossover of the conductance tensor in regions typically not accessible in the 
bosonic formalism. 


III. DUAL FERMIONIC VARIABLES AND RENORMALIZATION GROUP APPROACH TO THE 
CALCULATION OF THE CONDUCTANCE AT A JUNCTION OF THREE SPINFUL INTERACTING 

QUANTUM WIRES 

We now consider a three-spinful-wire junction, such as the one we sketch in Fig. [5] Resorting to the appropriate 
fermionic variables, we generalize to strongly interacting regions the weak-coupling FRG-approach. As a result, we 
map out the full dependence of the conductance tensor on the low-energy running cutoff scale even in strongly- 
interacting regions of the parameter space. Eventuall y, w e discuss the consistency of our results about the phase 
diagram of the junction with those obtained in Ref. [23j, particularly showing how our technique can be used to 
recover informations that typically cannot be derived within the bosonization approach used there. Consistently 
with Ref. [23j], in the following, we make the simplifying assumption that, in the weakly interacting regime, the bulk 
interaction is purely intra-wire and is the same in all the three wires. In fact, while this assumption is already expected 
to yields quite a rich phase diagram^, in principle our approach can be readily generalized to cases of different bulk 
interactions in different wires, such as the one considered in Ref. [&!]. 


A. The weakly interacting regime 


For a weak boundary interaction, assuming total spin conservation at the junction, the most relevant boundary 
interaction Hamiltonian is a combination of the bilinear operators in Eq. 0 . The relevant scattering processes at 
the junction are all encoded in the single-particle S'-matrix elements, ^ Because of spin conservation, 

the S'-matrix is diagonal in the spin index, that is, ,<r'){k) = fia,cr'Sjj'(k). Assuming also that the boundary 

Hamiltonian is symmetric under exchanging the wires with each other, the 3x3 matrix S(k) takes the form 


S = 


r t t 
t r t 
t t r 


(30) 


(Note that, in Eq. (1301) we assume that all the amplitudes are computed at the Fermi level and accordingly drop the 
index k from the S-matrix elements. This is consistent with the discussion of Appendix [A] where we assume that, 
close to the Fermi level, the scattering amplitudes are smooth functions of k.) It is worth mentioning that, in writing 
Eq. (1301) . we allowed for time-reversal symmetry breaking as a consequence, for instance, of a magnetic flux (j> piercing 
the centre of the junction (see Fig. 0. This implies that, in general, the scattering amplitude t from wire j to wire 
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j + 1 is different from the one from wire j to wire j — 1 (t). From Eq. (1301) one therefore finds that the F-matrix 
elements in Eq. (IA8I) are given by Fjji = § rSjj', with again /? = [~Si|| — 9i± + <? 2 ||] • On applying the FRG 

formalism of Appendix [Aj one readily obtains the RG-equations for the independent iS-matrix elements, given by 


dr 

dl 

dt 

dl 

dt 

dl 



\r\ 2 r — 2 ttr* 



* 




( 31 ) 


which, again, must be supplemented with the RG-equations for the running coupling strengths, Eqs. (1A10I11D . Since, 
as a consequence of spin conservation in scattering processes at the junction, the S'-matrix is diagonal in the spin 
indices, the charge- and spin-conductance tensors are equal to each other at the fixed points, as well as along the 
RG-trajectories obtained integrating Eqs. (13T1) . In particular, using the formalism of Appendix 1C 11 one obtains 


G c ,s(D) = — 

7r 


-R(D) + 1 
-T{D) 

- T(D) 


-T(D) 
—R{D) + 1 
- T(D) 


~T(D) 

- T(D) 
—R(D) + 1 


(32) 


with R(D) = |r(I?)| 2 , T(D) = |t(F)| 2 , and T(D) = \t(D)\ 2 . The RG-flow of the scattering coefficients T(D),T(D) is 
recovered by solving the set of differential equations 

^=-|[(5r-T)(l-T-T)-IT] 

^=-f[(5 f-T)(l-T-f)-TT\ , (33) 

which are derived from Eqs. (1311) by taking into account the unitarity constraint T(D) + T(D) + R(D) = 1. The fixed 
points of the boundary phase diagram are, therefore, determined by setting to zero the terms at the right-hand side 
of Eqs. (i33l) . From Eq. (l32l) one may therefore recover the corresponding charge- and the spin-conductance tensors. 
Borrowing the labels used in Ref. [23|, we obtain the following fixed points: 

• The[N c ,N s ] (’’disconnected”) fixed point 

This fixed point corresponds to having R = 1 and T = T = 0 which, according to Eq. (l32l) . yields 


2 0 0 0 

G C = G S = — 000 

n 0 0 0 


as it is appropriate for a disconnected junction. 

• The x++ fixed point 

This corresponds to R = T = 0, T = 1, which yields 


I Sw 


0 0 1 
10 0 
0 1 0 


and, accordingly 


G c = G 


s — 


7r 


10-1 
-1 1 0 
0-11 


(34) 


(35) 


(36) 


• The x _ fixed point 

This corresponds toi? = T = 0, T = 1, which yields 
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ln(D 0 /D) 


ln(D 0 /D) 


FIG. 6: Renormalization group flow of the scattering coefficients for a three-wire junction for gi t ± = 0, /3 constant and equal 
to -0.3, and different choices of the initial values of the scattering coefficients: 

(a) Renormalization group flow corresponding to T(Do) = T(Dq) = 0.1, R(Do) = 0.8. The curves corresponding to T(D) and 
to T(D) vs \n(Do/D) collapse onto the single red curve of the graph, while the flow of R{D) vs \n(Do/D) is described by the 
blue curve. As Do/D grows, the scattering coefficients flow towards the asymptotic values corresponding to the M-fixed point; 

(b) Renormalization group flow corresponding to T(Do) = 0.1,T(Do) = 0.3, R(Do) = 0.6. As Do/D grows, the scattering 

coefficients flow towards the asymptotic values corresponding to the x _-fixed point. 


and, accordingly 


I Sj,j 


0 1 0 
0 0 1 
10 0 



1-10 
0 1 -1 

-1 0 1 


(37) 


(38) 


• The M fixed point 

This corresponds to T = T , = |,f?=i and has to be identified with the symmetric)^, or with the Griffith^ 
fixed point of a junction of three interacting wires. One obtains 


and, accordingly 


r l 4 4 n 



L 9 9 9 J 


r 8 4 4 "1 



9 9 9 J 


(39) 


(40) 


X±± must clearly be identified with the ’’chiral” fixed points of Ref. [ 23| , where time-reversal symmetry breaking 
is maximum, both in the charge and in the spin sector. Along the RG-trajectories connecting two fixed points, the 
conductance flow is determined by Eq. (l32l) . The topology and the direction of the RG-trajectories depend on both 
/3 and on the bare values of the S'-matrix elements, ft scales with £ as determined by Eqs. (IA10I11I) . which makes it 
necessary to resort to a full numerical integration approach. A set of simplified situations can be realized, however, 
where keeps constant along RG-trajectories. For instance, if g\.±(Do) = 0, Eq. (fill) implies that ft is constant. 
In this case, from Eqs. (l33l) one readily sees that, if /? > 0, the boundary flow is towards the IVlV-fixed point. At 

variance, if f3 < 0 and T(Dq) > (<)T(D o), the boundary flow is towards the x++ (x_)-fixed point. As an example 

of possible RG-trajectories that may be realized in this specific case, in Fig. [6] we plot T{D) and T(D) for constant 
and negative, while we draw similar plots in Fig. [3 for constant and positive ft and in Fig. [5] for non-constant ft (see 
the captions for details). 

All the analysis we have done so far applies to a junction of three spinful quantum wires for weak bulk interaction. 
We now employ the dual-fermion approach to generalize the FRG-technique to regimes corresponding to strong bulk 
interactions either in the charge-, or in the spin-channel (or in both of them). 
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b) 


T(D)=T(D) 


R (D) 



ln(D 0 /D) 


ln(D 0 /D) 


FIG. 7: Renormalization group flow of the scattering coefficients for a three-wire junction for gi t ± = 0, /3 constant and equal 
to 0.3, and different choices of the initial values of the scattering coefficients: 

(a) Renormalization group flow corresponding to T(Do) = 0.35, T(Do) = 0.45, R(Do) = 0.2; 

(b) Renormalization group flow corresponding to T{Dq) = T(Dq) = 0.4, R(Do) = 0.2. As Do/D grows, in both cases the 
scattering coefficients flow towards the AW-fixed point. 




ln(D 0 /D) 


ln(D 0 /D) 


FIG. 8: Renormalization group flow of the scattering coefficients for a three-wire junction for non-constant /3: 

(a) Renormalization group flow corresponding to T(Do) = 0.2,T(Do) = 0.1,R(Do) = 0.4 and to <?i,|| (Do)/(2nv) = 
gi,±(Do)/(2TTv) = <? 2 ,|| {Do)/(2ttv) = gi t ±(Do)/(2nv) = 0.4. For these values of the bare parameters the junction is attracted 
by the disconnected fixed point (R —> 1 while T, T —> 0); 

(b) Renormalization group flow corresponding to T(Do) = 0A,T(Do) = 0.35, R(Do) = 0.25 and to gi t \\ (Do)/(2nv) = 
0.3, gi t ±(Do)/(2nv) = — 0.2, g 2 ,|| (Do)/(2nv) = gi t ±(Do)/(2nv) = 0.2. For these values of the bare parameters the junction 
is attracted by the x++-fi x ed point. 


B. Fermionic analysis of the strongly interacting regime at g c ~ g s ~ 3 

A first regime to which our dual-fermion approach can be successfully applied corresponds to a strong attractive 
interaction, both in the c harg e and in the spin channels. In particular, we assume g c ~ g s ~ 3. According to the phase 
diagram derived in Ref. [23( within the bosonization approach, in this range of values of the Luttinger parameters 
one expects to find a fixed point where paired electron tunneling and Andreev reflection are the dominant scattering 
processes at the junction and, in addition, two fixed points with maximally broken time-reversal symmetry, to be 

identified with the X++ an d the x_-fixed points discussed in the previous subsection. To apply the FRG-approach 

to this part of the phase diagram, we have to define the appropriate dual fermion coordinates. To do so, let us set 
gc — gs = 3. We therefore note that, though, in general, the charge- and spin-velocities u c and u s can be different 
from each other, one may easily make them equal by a pertinent rescaling of the real-space coordinate in the charge- 
and in the spin-sector of the bosonic Hamiltonian in Eq. (IB 101) . As the rescaling does not affect the boundary 
interaction (which is localized at x = 0), in the following, without loss of generality, we will assume u c = u s = u. In 
choosing the appropriate dual fermion coordinates, we use the criterion of mapping the fixed point we recover in the 
strongly interacting limit one-to-one onto those of the phase diagram in the weakly interacting regime. Referring to 
the bosonization formulas of Appendix [Bj we define the dual bosonic fields 4> c ( s )y (x), 0 C ( S ) j{x) in terms of those in 
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Eqs. (1B1IB2I) as 

’ fcW.lW ’ 

fc(*), 2 (a:) 
fc(*),3(®) 

®c(s),l ( X ) 

®c(s),2(-e) 

- &c(s),3( x ) - 

Consistently with Eqs. (El, we therefore define the dual fermionic fields as 

4>L,a,j(x) = 


111 

\/3 \/3 \/3 

\/3 \/3 \/3 

V3 \/3 75 

0 -1 


1 

-1 


0 

_ 1 
3 


1 

-1 

0 


1 _ 1 
3 3 

0 5 

s -I o 
1 1 1 
3\/3 3 >/3 3\/3 

?¥? 
3y/3 3V3 3V3 



' ^c( 5 ),lW ' 


$ c( s ),2(a;) 


$c(s),3(®) 


®c(s),l (*^) 


@c(«), 2 {x) 


- ®c(s), si.x') - 


(41) 


(42) 


Using Eq. (TUI) one sees that, when expressed in terms of the <f> and of the 0-helds, the bulk Hamiltonian in Eq. 
(IB10D reduces back to the one with g c = g s = 1, which, when expressed in terms of the fermionic fields defined in 
Eqs. CH, corresponds to the free Hamiltonian given by 


Hq,f = ~iu VV / dx \ ^ R , <r (x)d x 'ipR,j >e r(x) - $ L , a {x)d x i> LtjtCr {x) } . (43) 

- Jo 1 J 

Equation (14311 is the striking result of our technique of introducing dual fermion operators: it is a free-fermion 
Hamiltonian which describes a system that is strongly interacting in the original coordinates. Based upon the dual 
fermion fields in Eqs. (l42l) one may therefore introduce dual boundary operators analogous to those defined in Eq. 
namely, one may set 


0) = V’x j >(°)yx',j',cr(0) 


(44) 


and assume that the boundary interaction is realized as a linear combination of the operators in Eq. m and/or 
of products of two of them. In the absence of additional bulk interaction involving the dual fermion fields, or in 
the weakly interacting regime, the most relevant boundary interaction term is realized as a linear combination of the 
H-operators only. Therefore, the physically relevant processes at the junction are all encoded within the single-particle 
5-matrix elements in the basis of the dual fields, S(j.a)\(j' .a')- A nontrivial flow for the 5-matrix elements is induced 
by a nonzero bulk interaction in the dual-fermion theory, that is, by having g c ^ = 3 + 8 g c ^, with |(b?c(s) |/<? c (s) "C 1- 
The dual interaction Hamiltonian, H lnt can be readily recovered using Eqs. (1411421) . The result is 


■’ ('L ■’ pL 

flint = EE 9j;(cr,a') I &X PRj^{x)pL,j,o t'(%) E E / dxp R j,a(x)p.L,j',<r'(x) , (45) 

j=l CT,t7' 0 j^j '=1 er,(T' 0 


with pR( L ) tjt(T (x) =: i’R( L)>jt ' T (x)4>R(L),j,< t(x) and 

2iru(5g c + 6g s ) 2nu(8g c -8g s ) 

9j;(tT,<T') — g °cr,cr' + g 0<T,(T' 

87 Tu(Sg c + Sg s ) s 8ivu(6g c - 8g s ) s 

9(j,j');(cr,cr') — Q Q d cr,<?' > 


(46) 


plus terms that do not renormalize the scattering amplitudes. H lnt takes the form of the generalized bulk Hamiltonian 
in Eqs. (IA13IA14I) . Given the corresponding F-matrix elements reported in Eq. (1A16I) . one may derive the RG- 
equations in the case in which the spin is conserved at a scattering process at the junction, which implies <S'(j,< 7 );(j',<t / ) = 
da.a'Sj.j> , and the boundary interaction is symmetric under a cyclic permutation of the three wires, that is, the Sjj>- 
matrix elements are given by 


5 = 


fit 
t f t 
t t f 


(47) 
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Summing over both the inter-wire and the intra-wire processes allowed by the bulk interaction, one obtains 


df 

M 

dt 

dl 

dt 

d£ 


(7_a) _ |-|2- _ 2itr*} 

-{ 2 ^)mf\ 2 + ^ 2 n 

-(^) W |2 + W 2 r* } 


(48) 


with a = 5 , (jj') ; (cr,o-)/(27ru),7 = gj.^^/finu). Equations (l48l) are equal with Eqs. (1611) in the weakly interacting 
case, provided one substitutes the (running) parameter /3 in Eqs. (1311) with the (constant) parameter 7 — a. As a 
result, the RG-flow of the S-matrix elements is the same as the one obtained for the S'-matrix elements. Nevertheless, 
due to the nonlinear correspondence between the original and the dual fermionic fields, the result for the conductance 
at corresponding points of the phase diagram is completely different. To discuss this point, let us write the current 
operators at fixed spin polarization, J JiCr (a;), in terms of the dual fermionic fields as 


Jj,a{x)=eu ^2 {px,j-l,a(x) - px,j+l,a(x)} 

X=L,R 


(49) 


with j + 3 = j, and use the formalism of Appendix 1C 21 to derive the conductance tensor. From Eq. (IC19I) one 
eventually obtains 


G C (D) 


G S (D) = G 0 — 


3e 2 

^ T{D) 


(50) 


with 


'4 -2-2' 


C f(D)+f(D )) 

~T(D) 

- T(D ) 

-2 4 -2 

, r(D) = 

-T(D) 

C T(D)+T(D )) 

- T(D ) 

1 

tc 

1 

tc 


-T(D) 

~T{D) 

C AD) + T(D)) 


and T(D) = |t(D)| 2 , T(D) = \t(D)\ 2 . The RG-flow of T(D ) and of T(D) is determined by Eqs. (1531) . with /? replaced 
by 7 —a oc 5g c +5g s . As a result, for 5g c +5g s > 0 the stable fixed point is the ’’dual” disconnected fixed point, which we 
dub [N c , N s ], as, in bosonic coordinates, it corresponds to imposing Neumann boundary conditions on all the $ JjC ( s )(a;)- 
fields at x = 0. From Eqs. (150I51[) one therefore finds that the corresponding charge- and spin-conductance tensors 
are given by G c = G s = Gq. This is absolutely consistent with the result provided in Ref. [23| for g c = g s = 3. Indeed, 
from Eqs. (1411421) one sees that, resorting back to the original bosonic fields, the [N c> iV s ]-fixed points corresponds to 
the [D c , D s ]-fixed point of Ref. [2S| . with Dirichlct boundary conditions imposed on the relative fields <pi tC (s){x) = 
^[^i,c(s) (x) - $ 2 ,cO)(®)] and <p 2 ,c(s)(x) = ^g[$i, c (s)(z) + ^ 2 ,c(s)(*) - 2 $ 3 >c(s )(x)], where the charge- and the spin- 
conductance tensors for g c = g s = 3 are equal to each other and both equal to Go- To double-check the consistency 
between our dual-fermion FRG-formalism and the bosonization approach, we note that, at the [IV C , N s ]~ fixed point, a 
generic linear combination of the boundary operators in Eq. (1441) can be expressed, in the original bosonic degrees of 
freedom, as a linear combination of the operators Oy CT ( 0 ) = gLp,aPL,i.a e _ 5 [ e i+i,c(o)+ 03 , c (o)]- i f[ej + i, s (o)+ej, s ( 0 )] anc j 
of their Hermitean conjugates, which is the result obtained in Ref. [23f by means of a pertinent application of the 
delayed evaluation of boundary conditions (DEBC)-technique22^i. 

When 6 g c + Sg s < 0, the [N c , A^ s ]-fixed point becomes unstable. As in the weakly interacting case, we see that, if 

T(Dq) ^ T(Dq), the junction flows towards either one of the ”dual-chiral” fixed points, x++, x _, with the crossover 

of the conductance tensors with the scale being given by Eq. (l50l) . In particular, if T(D 0 ) > T(D 0 ), the flow is 
towards the infrared stable x ++ -fixed point. This corresponds to R = T = 0 ,T =1. From Eq. dSUl) . one therefore 
obtains that the fixed point conductances are given by 


G C = G 


S 


7 r 


1 -2 1 
1 1 -2 

-2 1 1 


e 2 

= —Q 

7 r 


(52) 


By consistency, one would expect that the X++ fixed point should be identified with the x++ fixed point emerging 
from the weak interaction calculation of the previous section. However, in order to compare the conductances obtained 
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in Eq. (l52l) with those of Eq. (l36l) one has to take into account that formula for the conductance tensor derived 
within dual-fermion approach applies to a junction connected to reservoirs with g c = g a = 3. Therefore, to make 
the comparison, one has to trade Eq. (1521) for a formula for the conductance tensors of a junction connected to 
reservoirs with g c = g s = 1, G c - w i,G s - tW i. As discussed in Appendix 1C 21 this can be done by using Eq. (IC16I) with 
Gin c(s) = ~ 9c<s \ —. The result is 

In ’ c W Sc(s)-1 w 


e 
7r 




i 

-i 

0 

1 

-1 ■ 
0 

l 3 / 

7r 

0 

-1 

1 


(53) 


that is, the same result as in Eq. (|36l) . Similarly, one can prove that the chiral %_-fixed point, towards which the 

RG-trajectories flow if T{Dq) < T(Dq), has to be identified with the %_-fixed point of Sec. IIII Al 

When T(D 0 ) = T(D 0 ), the RG-trajectories flow towards a nontrivial fixed point, which we dub M, by analogy to 
the M-fixed point we found in Sec. IIII Al Such a fixed point corresponds to R = 1/9, T = T = 4/9. Therefore, from 
Eq. (1501) . one finds that the fixed point conductance tensors are given by 


r 4 2 2 1 



L 3 3 3 J 

Performing the same transformation as in Eq. (USD, one eventually finds 



(54) 


(55) 


On comparing Eq. (1551) with Eq. (1401) we now see that, at odds with what happens with the chiral fixed points, the 
M and the M-fixed points cannot be identified with each other. While we are still lacking a clear explanation for this 
different behavior at different fixed points, we suspect that this shows that, while the conductance at X++ as well as the 
X —-fixed points are in a sense universal, that is, independent of the Luttinger parameters (provided one pertinently 
takes into account the corrections due to different Luttinger parameters for the reservoirs), the conductance at the 
M-fixed point does depend explicitly on the Luttinger parameters. This would definitely not be surprising, as such a 
feature would be shared by a similar fixed point such as, for instance, the nontrivial fixed point at a junction between 
a topological superconductor and two interacting one-dimensional electronic systems^. In any case, we believe that 
this issue calls for a deeper investigation, which will possibly be the subject of a forthcoming work. 

As so far we mainly concentrated around the ’’diagonal” in Luttinger parameter plane, that is, at g c ~ g s , we are 
now going to complement our analysis by discussing the regime with g c ~ 3,p s ~ 1, together with the complementary 
one, g c ~ l,g s ~ 3. 


C. Fermionic analysis of the strongly interacting regime for g c ~ 3, g s ~ 1 and g c ~ l,g s ~ 3 


We now discuss the ’’asymmetric” regime g c ~ 3,g s ~ 1. In order to recover the whole procedure, we again note 
that it is always possible to separately rescale the real-space coordinate in the bosonic Hamiltonian in Eq. (IB 101) . so 
to make the charge- and the spin-plasmon velocities to be both equal to u. Therefore, to actually define the dual 
fermion coordinates, let us assume g c = 3 ,g s = 1. Due to the absence of bulk interaction in the spin sector, we have 
no need to transform the Q s ,j- fields. At variance, we do trade the fields <d c j for the fields $ c ,j, @c,j defined 
in Eq. (TUI) . Accordingly, we consistently define the dual fermionic fields as 


XR,<r,j(x) = %,a,je5 [ ^’ c(x)+ ^> (x)+ ' T( ^' s(x)+0 ^' i(a:))] 

XL crj(x) = r}L J e5[* > 7’°^ _ ®7'c( a: )+ Cr ( $ 3,»( ;!: }“®3,»( x ))] 


(56) 


Again, one sees that, when expressed in terms of the dual fermion operators in Eqs. (1561) . the bulk Hamiltonian in 
Eq. (IB10I) reduces back to the free fermionic Hamiltonian, given by 
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Ho,f ;x = ~ iu 2^ 
j=i 
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L dx 


(57) 
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Having defined the dual fermion operators, we now assume that the leading boundary perturbation is realized as a 
linear combination of the dual boundary operators defined as 


A')(0) — Xxj,a(Q)XX',j',a(0) 


(58) 


and/or of products of two of them. Just as we have done before, we also assume that the most relevant scattering 
processes at the junction are fully described by means of the single-particle S'-matrix elements in the basis of the 
X-fields, Slightly displacing ( g c ,gs ) from (3,1), that is, setting g c = 3 + 5g c ,g s = 1 + 5g s , with 

|<5<7 c |/3, |<$<7s| "C 1; gives rise to an effective interaction Hamiltonian H x - i n t, which takes exactly the same form as H lnt 
in Eq. (l45l) . and is given by 


!int. — 




j = 1 <T,tr' 


/ d x Px\R jV ( x )Px; L 
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Px'iRiji* 7 Px’^ij'i G ' (*^) j 


with p x;R (L),jA x ) = : Xr( L )jA x ) Xi and 


_ 87 Tu(Sg c - 35g s ) x i 87 iu{5g c + 3 5g s ) s 

9 x\r,(a,°') ~ 9 + 9 °°,a' 

Anu(5g c -35g s ) x 4wu(5g c + 35g s ) r 

~ g dcr > CT ' g > 


(60) 


plus terms that do not renormalize the scattering amplitudes. Making the assumption that the spin is conserved at a 
scattering process at the junction, we again obtain that S x -,(j, 0 );(j' ,<?’) = '-boO'd')’ w ith, for a boundary interaction 

symmetric under a cyclic permutation of the three wires, the ^-matrix being given by 


Sx = 


T ^ tx 

he r X tx 
tx t-X r x 


(61) 


The RG-equations for the running S^-matrix elements are derived in perfect analogy with Eq. (1481) . The result is 
exactly the same, except that now 7 — a oc 5g c — 35g s . In order to trace out the correspondence between the fixed 
point of the boundary phase diagram for the dual-fermion scattering amplitudes and those of the phase diagram for 
the original fermion amplitudes, we now discuss the behavior of the charge- and of the spin-conductance tensor along 
the RG-trajectories. To do so, we note that, due to the fact that the spin sector of the theory is left unchanged, when 
resorting to the dual coordinates, the spin-conductance tensor, when expressed in terms of the scattering coefficients at 
the junction, takes the same form as in the noninteracting case, given in Eq. (1321) . As variance, the charge-conductance 
tensor depends on the scattering coefficients as given in Eq. ©. As a result, one obtains 


3e 2 

Gc(D) = Go - —r X (D) 


G e (D) = - 

7 r 


~Rx(D) +1 -T X (D) -T X {D) 

- TJD) —R\ (D) + 1 -TJD) 

~T X {D) - T X (D) -R X (D) + 1 


(62) 


with R X {D) = \r x (D)\ 2 ,T X (D) = \t x (D)\ 2 , T X (D) = \t x (D)\ 2 . We are, now, in the position of mapping out the whole 
phase diagram of the spinful junction for g c ~ 3 ,g s ~ 1 , including the fixed point manifold, and of tracing out the 
correspondence between the fixed points given in terms of the dual fermion amplitudes, and the described in terms of 
the original fermionic coordinates^. First of all, we note that, when 5g c — 3 6g s > 0, the system is attracted towards 
the ’’dual disconnected” fixed point [A ^ x , c ,N XtS ], characterized by the scattering coefficients R x = 1,T X = T x = 0. 
At such a fixed point, one obtains G c = Go,G s = 0, which enables us to identify [N XtC ,N x ^ s ] with the [D c , iV s ]-fixed 
point in the phase diagram of Ref. [23|, that is, with a spin-insulating fixed point where the most relevant process at 
the junction is pair-correlated Andreev reflection in each wire. At variance, when 5g c — 3 5g s < 0, the junction flows 

towards one among the dual x++,X _> or M- fixed points. In particular, from Eqs. (l62l) . one sees that, at the dual 

X_l_+-fixed point, G c is given by Eq. (l5!?l) . while G s takes the form provided in Eq. (15S1) . After the correction of Eq. 
ED , one eventually finds that G c and G s are equal to each other, and both equal to the fixed-point conductance 
at the x ++ -fixed point in the original coordinates. Thus, we are eventually led to identify the dual x + _|_-fixed point 
with the analogous one, realized in the original coordinates. A similar argument leads to the identification of the dual 

X _-fixed point with the analogous one, realized in the original coordinates. As for what concerns the dual M-fixed 

point, after correcting the charge-conductance tensor as in Eq. (l55l) . one finds that, at such a fixed point, 
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Putting together Eqs. (I63I55I40I) we again see that M-like fixed points are not mapped onto each other, not even 
after the correction of Eq. (1551) . This is again consistent with our previous hypothesis, namely, that the conductance 
at the M-fixed point does depend explicitly on the Luttinger parameters and, therefore, it is different in the various 
cases we discussed before. 

Before concluding this subsection, we point out that the same analysis we just performed in the case g c ~ 3, g s ~ 1 
does apply equally well to the complementary situation g c ~ 1, <? s ~ 3, provided one swaps charge- and spin-operators 
(and conductances) with each other. 

Putting together all the results we obtained using the FRG-approach, one may infer the global topology of the 
phase diagram of a spinful three-wire junction and compare the results with those obtained within the bosonization 
approach. This will be the subject of the next subsection. 


D. Global topology of the phase diagram from fermionic renormalization group approach 

Following Ref. [23J, we discuss the main features of the phase diagram of the three-wire junction within various 
regions in the g c — g s plane. Let us start from the ”quasisymmetric” region g c ~ g s . From the results of Sec. IIII Al 
we see that, setting g c = 1 + Sg c ,g a = 1 + Sg s , as long as Sg c + Sg s < 0 (corresponding to + -A- < 1), the system 

flows towards the disconnected [N c , 7V s ]-fixed point. At variance, for Sg c + Sg s > 0 (that is, for > 1), as 

soon as scattering processes from wire j to wires j ± 1 take place at different rates (i.e., T ^T), the RG-trajectories 

flow towards either one of the X++ or X_-fixed points. While this is basically consistent with the region of the 

phase diagram derived in^ corresponding to g c ~ g s ~ 1, in addition, when T = T, we found that the system flows 
towards a symmetric fixed point, which we dubbed iff, with peculiar, g c , g s -dependent transport properties. Within 
the FRG-approach we were able to map out the full crossover of the charge- and spin-conductance tensor between 
any two of the fixed points listed above, with some paradigmatic examples shown in the figures of Sec. IIII Al Keeping 
within the quasisymmetric region, in Sec. IIII Bl we show that, setting g c = 3 + Sg c ,g s = 3 + 5g s , for Sq r + Sg s > 0 
(that is, for g c + g s > 6), the stable RG-fixed point corresponds to the [D c , D s ]-fixed point of Ref. |23|] while, as 
soon as 5g c + Sg s < 0 (that is, for g c + g s < 6), the system flows towards either one of the x++ or X— fixed points 
in the non-symmetric case, or towards an ” M-like” fixed point in the symmetric case. As discussed above, while, 

at both the x++ and the x _-fixed points the conductance tensors for the junction not connected to the leads are 

the same, regardless of the value of the Luttinger parameters, at variance, at the M-fixed point they do depend on 
g c and g s and, in this sense, they appear to be ” nonuniversal”. Over all, the results we obtained across the region 
9 c ~ gs are consistent with a phase diagram where the [N C ,N S \ and the [D C ,D S \-fixed points are respectively stable 
for A—f 2 g < 1 and for g c + g s > 6, while, at intermediate values of the Luttinger parameters, depending on the 

bare values of the scattering coefficients at the junction, one out of the (universal) chiral X++ • X_-fixed points or 

the (nonuniversal) M-fixed point becomes stable. This results already complements the phase diagram of Ref. [23j 
by introducing the M-fixed point, which has necessarily to be there, in order to separate the phases corresponding to 

X++ and to x _from each other. To push our analysis outside of the g c ~ g s -region, we discussed the nonsymmetric 

case g c ~ 3 ,g s ~ 1. In this case, we found that the manifold of fixed points consists of the [D C ,N S ] asymmetric 
fixed point, at which the junction is characterized by perfect pair-correlated Andreev reflection in each wire, while 

it is perfectly insulating in the spin sector, the chiral x++,X _-fixed points and, again, an M-like fixed point. Also 

in this region our results appear on one hand to be consistent with the phase diagram of Ref. [23j . on the other 
hand to complement it with singling out the M-like fixed point. The complementary regime g c ~ l,g s ~ 3 can be 
straightforwardly recovered from the previous discussion by just swapping charge and spin with each other. Aside 
from recovering the global phase diagram of the junction, our technique allows for generalizing to strongly-interacting 
regimes the main advantage of using fermionic, rather than bosonic coordinates, that is, the possibility mapping out 
the crossover of the conductance between fixed points in the phase diagram. 


IV. DISCUSSION AND CONCLUSIONS 

In the paper, we generalize the RG approach to junctions of strongly-interacting QWs. In order to do so, we make 
a combined use of both the fermionic and the bosonic approaches to interacting electronic systems in one dimension, 
which enables us to build pertinent nonlocal transformation between the original fermion fields and dual-fermion 
operators, so that, a theory that is strongly interacting in terms of the former ones, maps onto a weakly interacting 
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one, in terms of the latter ones. On combining the dual-fermion approach with the FRG-technique, we are able to 
produce new and interesting results, already for the well-known two wire junction. When applied to a Z 3 -symmetric 
three-wire junction, our technique first of all allows for recovering fundamental informations concerning the topology 
of the global phase diagram, as well as all the fixed points accessible to the junction in various regions of the parameter 
space. While, in this respect, our approach looks like a useful means to complement the bosonization approach to 
conductance properties of junctions of quantum wires, where it appears extremely useful and, in a sense, rather 
unique, is in providing the full crossover of the conductance tensors between any two fixed points connected by an 
RG-trajectory. The crossover in the conductance properties can be experimentally mapped out by monitoring the 
transport properties of the junction as a function of a running reference scale, such as the temperature, or the effective 
system size. While the standard FRG approach just yields crossover curves at weak bulk interaction in the quantum 
wirea^ ~ 15 i 18 , as stated above, our approach extends such a virtue of the fermionic approach to regions at strong values 
of the bulk electronic interaction in the wires. 

By resorting to the appropriate dual fermionic degrees of freedom, our approach allows for describing in terms of 
effectively one-particle S'-matrix elements correlated pair scattering and/or Andreev reflection, in regions of values 
of the interaction parameters where they correspond to the most relevant scattering processes at the junction. This 
allows for envisaging, within our technique, fixed points such as the [D c , D s ], or the [D c , iV s ] one. In fact, due to basic 
assumption of the standard FRG-approach that all the relevant processes at the junction are encoded in the single¬ 
particle S'-matrix elements, fixed points such as those listed before are typically not expected to be recovered without 
resorting to the appropriate dual fermion coordinates, not even after relaxing the weak bulk interaction constraint^. 

While, for simplicity, here we restrict ourselves to the case of a symmetric junction and of a spin-conserving boundary 
interaction at the junction, our approach can be readily generalized to a non-symmetric junction characterized, for 
instance, by different Luttinger parameters in different wires^*, and/or by a non-spin-conserving boundary interaction. 
Also, a generalization of our approach to a junction involving ordinary^, or topological superconductor s 71 ' 72 is likely to 
allow for describing the full crossover of the conductance in a single junction, as well as of the equilibrium (Josephson) 
current in a SNS-junction thus generalizing, in this latter case, the results obtained Refs. [73H75| to an SNS-junction 
with an interacting central region. 

Finally, it is worth stressing that our technique generically complements the RG-approach, so to extend it to 
strongly-interacting problems. Very likely, rather than combining it with the FRG-technique, one could work out 
the RG-trajectories by using the alternative (and, to some extent, more accurate, though less tractable analytically) 
fRG-approach^—. Although we believe this is a potentially interesting research topic to pursue, it goes beyond 
the scope of this work, where, for the sake of simplicity and analytical tractability, we rather preferred to use the 
FRG-technique. 

We would like to thank I. Affleck and Z. Shi for enlightening discussions and for sharing with us their unpublished 
results about the FRG-approach to a junction of three spinless quantum wires. We acknowledge insightful discussions 
with A. Zazunov and E. Eriksson at various stages of completion of this work. A. N. acknowledges financial support 
from European Commission, European Social Fund and Regione Calabria. 


Appendix A: Derivation of the fermionic renormalization group equations for the S'-matrix 


In this Appendix we review the derivation of the FRG-equations for the S'-matrix elements describing single-particle 
scattering at a junction of quantum wires, as discussed in Refs. mm- In our paper we compute dc-transport 
properties of junctions of quantum wires. In doing so, we describe each wire by only retaining its low-energy, long- 
wavelength excitations about the Fermi points ztkp. Finally, we alledge for the system to present an ’’inner” boundary 
at x = 0, where the junction is located and the boundary conditions on the fields are determined by the boundary 
interaction describing the junction, and an ’’outer” boundary at x = L which is just required for introducing a cutoff 
length scale L, eventually sent to oo at the end of the calculations. As a result, the ’’bulk” Hamiltonian for the wires 
is realized as Ffsuik = Ho + Hy, with Hq being the free Hamiltonian for noninteracting chiral fermions in Eq. Q, 
while the two-body interaction Hamiltonian Hy is given by 



L 


dx dy pj^(x)Vj. a ^{x 


y)pj,*' (y) 


(Al) 


Note that, as typically done in this class of problems^ 15 d 8 ; i n Eq. im we are assuming a purely ”intra-wire” 
interaction. In fact, in the analysis of our paper, we had to deal with an emerging bulk interaction Hamiltonian with 
nontrivial inter-wire interaction terms. Yet, for the sake of simplicity, we prefer to add inter-wire interactions to the 
simplified model Hamiltonian, where a local approximation for the interaction potential is done. Indeed, on assuming 
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a short-range interaction potential effective over a typical length scale A (as it happens, for instance, in gated wires), 
one may approximate Eq. ED by setting x ~ y in the Coulomb potential, so that Hy is approximated by H lnt in Eq. 
([2]) . In order to derive Eq. (j2|) . we have neglected terms obtained integrating operators proportional to the rapidly 
oscillating functions e ±2lfcFa: . Based on an analogous argument, we assume that the umklapp term (oc e ±4lkFX ) does 
not appear either (in fact, an additional argument, related to the irrelevance of the corresponding operator, can be 
recovered within the bosonization approach to the proble m 76 ’ 77 ). Finally, it is worth mentioning that we also ignore 
terms of the form 


Hint,g 4 — 2 E/ E/ Vj,acr'{ 0) 

3 =1 a a' 

since they just renormalize the Fermi velocity and the chemical potential^ and, so, do not effectively contribute to 
the renormalization of the ^-matrix at the junction. Two relevant limiting cases can be recovered from Eq. (|2]l. The 
former one corresponds to spinless fermions and is recovered by dropping terms containing fields with a given spin 
polarization, say a =j.. In this case, on setting t it easy to verify that Eq. © reduces to 


dx [ pL,j,a(x)pL,j,cr'(x ) + PR,j,a{x)PR,j,<j' (x)] 


(A2) 
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(A 3 ) 


in agreement with Eq. (6) of Ref. [lsj . The latter one corresponds to assuming spinful fermions, but a spin independent 
interaction, setting < 7 yi( 2 )|| = ffy 1 ( 2 ),j_ = 9j, 1 ( 2 )- In this case, Eq. © reduces to 


Hn 
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dx 


9jA^RJ,A x )i’L,j,a'( X )^RJ,^(x)ilLjA x )+9j^Rj,a(x)i’L, j ,a'(x)^ L ,j,tT'{x)i>R,j,a{x) , ( A4 ) 


in agreement with Eq. (50) of Ref. [l8j]- As already stated in the main text, FRG approach typically applies to the 
case in which the boundary scattering processes are fully described by the single-particle S'-matrix elements in Eq. 
Q. The first FRG step consists in trading the quartic bulk interaction Hamiltonian for a quadratic one, by means 
of a pertinent Hartree-Fock (HF) decomposition of the interactio n 13 ’ 14 i 18 . After performing the HF decomposition 
and making the standard assumption that, as we are dealing with low-energy excitations around the Fermi level, the 
S-matrix elements can be taken to be all independent of energy, H ln t reduces to 


jj, ^ fa, 2 , || 9j, i,|| 9j, i,-l) j L dx 

3 = 1 
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for spinful electrons, and to 
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(A6) 


for spinless electrons. The corrections to the amplitude for an incoming/outgoing electron in wire j to become an 
outgoing/incoming electron in wire j' under the effect of H m t are readily computed using the approach o f 13 ’ 14 ’ 18 . The 
corresponding corrections to the Sjji matrix elements are given by 
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dSw = - E (> (^ t ) i ,y Si'j. - d£ 


(A 7 ) 


with dl = dln(i/A), A being a length scale corresponding to the (finite) range of the interaction potential Vy CTCT '(x) 
and the Friedel matrix F defined as the block-diagonal matrix 
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Fjj' ~ ^ nv (ft,2,|| ft,i,|| ft,i,J-) Sjj6jj> 

in the spinful case, while 


(A8) 


Fj ’ j ' ~ 4 ttv (// '' 2 9jX ^ 


(A9) 


in the spinless case. The term ln(L/A) is due the presence of the x~ x in the integral giving the correction to the 
amplitudes within the HF-approximation. Throughout the poor man’s scaling method, one can replace dl = dln(L/A) 
with dln(iAo/-D), with Dq being a high-energy cutoff (~ bandwidth of bulk electrons) and D being the low-energy 
running scale. The RG-equations for the iS-matrix must be supplemented with the RG-equations for the interaction 
constants gi m, gi,±, g 2 ,\\, g 2 ,±- The procedure is discussed in detail, for instance, in Ref. [78[. Here, we just provide 
the final result for the scaling equations of the bulk interaction strengths, which isZ£ 


d 9j, 2,|| 
dl 

dgj, 2,_L 

dl 

d 9j. i,|| 
dl 

dftp ,± 

dl 

valid in the weak coupling regime. For a 


1 

2 irv 

1 

2ttv 


(ft, i,|| ) 2 

(ft, i,_l) 2 


- [(ftM.ll) 2 + (ft,F- l) 2 ] 

= _ ^ 2 ft,b-L [ft,2,1. - ft,2,|| + ft,l,||] » 

spin-independent interaction, Eqs. (1A10I) reduce to 
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with S = 0 for spinless electrons and S' = 1/2 otherwis e^ 13 ! 14 . Clearly, the renormalization group equations for the 
S-matrix must be supplemented with Eqs. (IA10I) in the spinful case in which, differently from what happens in the 
spinless case, there is a nontrivial flow of the bulk interaction parameters with the scale L/ A. 

In the analysis of our paper, we also have to consider a generalized bulk interaction which, for gj,i,x = 0, has a 
nonzero inter-wire component, that is, the generalized bulk interaction Hamiltonian reads 
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with 
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which reduces to Flint in Eq. 


Mj'=i 

provided one identifies the intra-wire interaction strengths as 


ft;(t,t) = ft;(Id) = ft,2.|| - ft,i,|| 

ft ;(t,T) = ft ;(T,t) = ft, 2,4. , (A15) 

and sets g(j,j')-la,n') = 0 for j ^ j'. On assuming again that the E-matrix takes the spin-diagonal form in Eq. ([7|i. one 
finds that the RG-equations for the E-matrix elements are again those provided in Eq. (lATl) . but with the E-matrix 
elements now given by 
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As we discuss in the main text, the modification in Eq. (IA16I) does not substantially affect the solutions of the RG 
equations. 
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Appendix B: Review of bosonization rules for interacting one-dimensional quantum wires 


In this Appendix, we review the basic bosonization formulas for one-dimensional interacting electronic systems, 
which provide us with the fundamental formal ground, on which we relate when defining the effective fermionic 
coordinates in the strongly interacting regimes. The first step towards a fully bosonic description of the junction is 
to rewrite in bosonic coordinates the free Hamiltonian Hq in Eq. ©■ This is done by introducing I\ bosonic fields 
4> c j, described by the (charge-sector) Hamiltonian 


K 


H, 


0 ;B;c 


-c/. *£ 


-(d t $ cJ ) 2 +v(d x $ CtJ ) 2 


3= 1 L 


and K bosonic fields described by the (spin-sector) Hamiltonian 


K 


Hq;B;s — 


47r 


dx 

i=i 


-(d t $ s ,j) 2 +v(d x $ s j) 2 


(Bl) 


(B2) 


together with the corresponding dual fields 0 c ,j, 0 s y, which are related to the <f>-fields by means of the cross derivative 
relations 


dtQj, c( s ) = vd x $ jtC(s ) . (B3) 

Therefore, one rewrites the chiral fermionic fields in terms of the bosonic coordinates introduced above as 

3 pR,a,j(x) = e 3 . (*))] 

il/L^iix) = T} L , a je* [ **' 4d-e ilC ( I )+dh>W-0i,»( I ))] , (B4) 

with rjR uj and rjL,o,j being real-fermion Klein factors. Note that the vertex operators in Eq. (IB4I) have been set 
consistently with the fact that, as a consequence of Eqs. (IB3I) . one finds that the following linear combinations are 
chiral fields 


4*R,c(s),j (•£> 0 “1“ ®j,c(s) (*£? 0] 

4 ) L,c(s),j (*^5 [ ^J,c(s) (^5 “h Oji,c(s) 0] 

In bosonic coordinates, the chiral fermionic density operators at fixed spin polarization are realized as 

: ' t l J R,a,j( x )' l pR,<r,j( x ) '■ = ^{-9x0j,c(a;) + 9a;$j iC (a;) + cr[-9 x 0j !S (a;) + S x $ J)S (a;)]} 

: '0Z, <r Ax)ipL,a,j(x) : = --^-{d x Q jtC (x) + d x $j. c (x) + a[d x e jtS (x) + 5 x $ JiS (a;)]} , 

’ o7T 


(B5) 


(B6) 


with the double columns :: denoting normal-ordering with respect to the fermionic groundstate. From Eqs. m one 
therefore recovers the charge- and the spin-density operators given by 


Pc,j{x) = y]{: i’R, er j{x)i>R,aj{x)+^l ><Tij {x)'ip Ltl7tj (x) :} —t d x e jtC {x) 

a 

Ps,j(x) = y]cr{: ^ aJ {x)i>R^j{x) + ^ L aj {x)^ Ly<JJ {x) :} d x O jtS {x) 

£7 

together with the corresponding current operators realized as 

Jc,j{x) = y]{: i’R,„,j(x)ll)R,< r ,j(x) - •tpL,*,j( x )' t l’ L ,<r,j(x) ■} ^ ^ &X&j,c( X ) 


(B7) 


Js,j(x) = ^R^jW^R^Ax) - A L ,aA X ^ L ^A X ) '■} 7,7 dx®j,s( x ) 


(B8) 
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All the previous transformations lead to contributions to the system Hamiltonian that are quadratic in the bosonic 
fields. The interaction Hamiltonian shares this property, once one has set to zero the all the terms oc in which 

case, on rewriting iAnt in Eq. (0 in terms of the bosonic fields, one obtains H lnt = H\ nt i + iAnt, 2 , with 


tt _ Si,i.ll gj,2,|| [ L 

-Clint,1 — 2 _^ 


j=1 

K 


167T 2 


dx {-(d x ®j, c {x)) 2 - (d x Oj tS (x)) 2 + (d x $ jtC (x)) 2 + (d x $j, s (x)) 2 } 


-Hint,2 — T. 


j=1 


gj,2,-L 

167T 2 


dx {(d x e jiC (x)) 2 + {d x $ jtS (x)) 2 - (d x Qj tS (x)) 2 - (d x $ jyC {x)) 2 } . 


(B9) 


On adding the terms in Eqs. (IB9I) to the noninteracting Hamiltonian for the charge sector (Eq. (IBljl 'l plus the one 
for the spin sector (Eq. (IB2l) h one eventually obtains a bosonic Hamiltonian H that is still quadratic, although with 
pertinently renormalized coefficients, given by 


1 ' / f 1 

H = — y Uj.c / dx < g jtC (d x $ jtC (x)) 2 4- (d x Q j c (x)f 

j~i l 9j,c 

1 f L ( 1 

+ “7 y ^ I dx \ (*^)) 

4?r J o l 9j,s 


(BIO) 


with 


Uj,c9j,c — ^ 
u. 




= V 


9j,c 

Uj,s9j,s = v 
Uj 




= V 


9j , 


1 - 


1 - 


gj,1,11 ~ gj,2,11 - gj,2,_L 
47 TV 

9j, 1,11 - gj,2,|| - gj,2,_L 
4irv 

gj,1,11 ~ flj.2,|| +gj,2._L 
47ft; 

gj,1,11 g j ,2,|| +gj,2,-L 

Aitv 


(Bll) 


(Note the use of Eqs. (IB3I) to express H in terms of both the $- and the O-fields.) When gyi.j_ ^ 0, one can employ 
the identities 


lj L,<r,j( x )i > L,e,3 

to express the total additional contribution to H lnt , H[ nt 3 , as 

,L K 

Hint ,3 ~ / dx y Wj cos[2$ iiS (x)] , 

Jo j=1 

with Wj oc gj.i,±- To lowest order, the renormalization group equation for Wj is 




4 ai (^L,aj(x) -O 


dW 3 

di 


2 - 


gj,S- 


Wf : 


(B12) 


(B13) 


(B14) 


which tells us that the nonlinear interaction term oc gj,i,± is irrelevant as long as gj }S < 1 and, accordingly, even 
if gj,i,-L has not been fine-tuned to 0, it can be dropped out of the effective low-energy, long-wavelength theory. In 
general, when using the bosonization, we assume that either this is the case, or that all the gj, i,j_ have been fine- 
tuned to 0. As for what concerns the boundary Hamiltonian Hb describing the junction between the quantum wires, 
throughout all the paper we have assumed that total spin was conserved at each scattering event at the junction. This 
basically implies that Hb can be fully as in Eq. © of the main text. While for a weak bulk interaction within the 
wires the most relevant contribution to Hb is typically realized as a linear combinations of (some of) the operators 
in Eq. ©, at strong enough attraction either in the charge-, or in the spin-channel (or both), as we discuss in the 
paper, the most relevant boundary interaction term can be realized as a product of two H-operators. In any case, 
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all the boundary operators can be rewritten in bosonic coordinates by pertinently employing the bosonization rules 
in Eqs. m- As the junction naturally defines a (common) boundary for all the quantum wires, Eqs. m must 
typically be supplemented with the appropriate boundary conditions for the <I> ? - and for the ©./-fields. This ’’delayed 
evaluation of the boundary conditions” (DEBC)-procedure can be typically implemented in the correspondence of a 
conformal fixed poin t 20 ' 21 . It has the net effect of making different U-operators in Eq. a to ’’collapse” onto each 
other, thus substantially reducing the number of independent boundary operators that can potentially enter Hb and, 
thus, strongly constraining the number of relevant scattering processes that can take place in the proximity of a given 
fixed point. For instance, the fields ^j,c(«)0c) corresponding to a wire that is disconnected from the junction must 
obey Neumann boundary conditions such as d :l Aj,c(s) (x = 0) = 0 and, correspondingly, the fields ©j, c (.s) (x) must 
be pinned at x = 0 (Dirichlet boundary conditions). Such boundary conditions only allow for nontrivial boundary 
operators in the form ^(j,j'), CT ,(j?,L)(0) ~ exp [—|($ c j(0) — d> c j/(0)) — ^($ 5J (0) — (0))] (plus the corresponding 

Hermitean conjugates). Various alternative possible situations are discussed in the paper. 


Appendix C: Linear response theory approach to the conductance tensor of a junction of quantum wires 

In this section, we review the linear response theory approach to the derivation of the conductance tensor for a 
junction of quantum wires. After developing the main formula for computing the conductance tensor within linear 
response theory, we will apply it to a number specific cases we analyze in the main text of our paper. 

Let Wj be the j-th wire connected to the junction (j = 1,..., K). In order to implement linear response theory, 
we imagine to connect Wj to a reservoir IZj , which can either be characterized by the same parameters as the wire 
to which it is connected, or not (a typical situation corresponds to Wj being an interacting one-dimensional quantum 
wire, described as a single Luttinger liquid, connected to a Fermi liquid reservoir—). For the sake of simplicity, we 
assume that the parameters characterizing both the wires and the reservoirs are all independent of j. In order to 
induce a current flow across the junction, we assume that each reservoir is characterized by an equilibrium distribution 
for particles with spin a , with chemical potential fJij,„ = eVj i<7 . This induces electric fields {£y iCr (i)} distributed in 
the various branches of the junction, which we account for by introducing a set of uniform vector potentials Aj. a (t), 
one for each wire, such that £j,a(t) = —d t A ha (t). Letting Jj tCr (x) be the current operator for particles with spin a 
in wire-j and letting each wire to be of length L , we may define the ’’source” Hamiltonian Hs OUIce {t ), describing the 
coupling to the applied electric fields, given by 


K 


H So 


■<*> = 

j= 1 a 


dx Aj,cr{t)Jj^{x) 


(Cl) 


Let us, now, denote with J^j a yi{x,t), Hs OUIce; j(t) the operators taken in the interaction representation with respect 
to the Hamiltonian without the source term. Within linear response theory, the current of particles with spin a 
evaluated at point x of wire-j is therefore given by 

K rOO pL 

I jA x A) = i X] / dt ' / dx ' V U,'ry,U','r')( X ’ t: ’ X 'A)Aj>,a'{£) , (C2) 

j'=l a’ 

with 


,a') (*£? ) ^(j' t ^)? J{j' )]) 


(C3) 


Equation (IC2I) does generically apply to any situation, whether the wires are interacting, or not, and whether one 
uses a fermionic, or a bosonic representation for the Hamiltonian of the junction. An important remark, however, is 
that, in any case, the current must be consistently probed outside of the region across which the electric fields are 
applied, that is, in Eq. (IC2I) one has always to assume that x > L. Another important formula is the Fourier-space 
counterpart of Eq. (IC2I) , that is 


Ij.a (■£; ^u) 



dx' V Uta) . (jV /) (x, a -!\u)£j’y (w) , 


(C4) 


with 
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= j 

f dt e lut Ij tCr (x,t) 

A 

|T 

ii 

f dt e lut £jy a '(t) 

ii 

'a 

[ dt e lut 'D {j ^ j y a , ) {x, t; x', 0) 


(C5) 


In the following part of this Appendix, we consider some specific applications to various cases we analyze in the 
main text. An important observation to make at this point is that, on connecting the wires to the reservoirs, we 
basically assume a continuity condition for the current operator at the interface. From the microscopical point of 
view, this appears to be the ’’macroscopic” counterpart of the ’’smooth” crossover in the interaction strength in real 
space discussed in the microscopic lattice model considered in Ref. j53j. It would be also interesting to work out 
a macroscopic field-theoretical Hamiltonian describing a ’’sharp” interface in the microscopic model, but this goes 
beyond of the scope of this work. As specified above, in the following we assume current continuity at the interface, 
that is, a smooth crossover in the bulk interaction strength. 


1. Junction of noninteracting quantum wires 


Due to the absence of multi-particle processes, the conductance tensor for a junction of noninteracting quantum wires 
can be fully expressed in terms of only single-fermion S-matrix elements. The bulk Hamiltonian for a noninteracting 
junction of K quantum wires is given in Eq. ©■ Within each wire the chiral fields can be expressed in terms of their 
Fourier modes as 

i> rjA x ) = etkXa RjA k ) 

4>L,jA X ) e ~ lkX(1 L,j,a{k) , (C6) 

with the right-handed and the left-handed chiral modes related to each other by the S'-matrix elements as 

K 

a R,j,<r(k) = ^ ( 'y ' {k) a L,j'{k) . (C7) 

j'= 1 <*' 

The current operator in wire-j for particles with spin polarization a is given by 

JjA x ) = ev i : ^R,jA x AR,iA x ) ■ - ■ ^l,jA x ^ l ^A x ) : l > ( C8 ) 

which, taking into account that the S-matrix is diagonal in the spin indices, implies 

')( x > x 'iA = e 2 v 2 6 <J ,cr'{'D(R,Ry,j,j' {x, x'\u)) + 'D(L,Ly,j,j' 0 , x'-,u>)- T>(R, L yj,j'{x, x'\w)- x'\u>)} . 

(C9) 

Assuming a thermal distribution for fermions in the reservoirs, one therefore obtains 


‘ D (R,Ry,j,j'(. x ’ x ']u) = 


V(L,Ly,j,j’(x,x'-,w) = 


^j,j' 

4n 2 v 2 

4n 2 v 2 


dE dE' f (E)f(E') 


dE dE’ f (E)f(E') 
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— i(^ E \ E x+x') i ^ E \ E ^ {x+x') 


Vf 


*1 s i'. 
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with Sj.ji denoting the S'-matrix elements computed at the Fermi level, f{E) being the Fermi distribution function, 
and 77 = 0 + . Since, as stated above, we assume that the current is measured outside of the region over which the 
electric fields are applied, in the following we assume x, x' > 0, 0 < x' < L and x > L. As a result, integrating in 
Eqs. (1C 101) over dE' and closing the integration path over the appropriate half-plane according to what is the phase 
factor in the argument of the integrals, one obtains 

J dEe^^-^f(E){f(-E-co)-fi-E + uj)} 

V {RM , h y{x,x'- ) u) = ^^- j dEeW x+x ')f(E){f(-E-w)-f{-E + «j)} 

/ D( L ,R)-,j,j'(x,x , ;u) = 0 . (Cll) 

Going to the zero -T limit and focusing onto the dc current (w —> 0), one recovers that the current in wire-j for 
particles with spin polarization tr, is given by limit 


j'=1 

The dc-conductance matrix elements , ff ') are therefore given by 

r)T p 2 

{ —|<Sj,j '| 2 + hi'} , 

The charge- and the spin-conductance tensors G C ,G S , are respectively defined as 


G, 


Mb') 




U,V);U',&') > G S ;(jJ') 


= V aa'G 


a,<j' 


(leOMW') ) 


(C12) 


(C13) 


(C14) 


which, from Eq. (1C13I) . implies 




(CIS) 


When a weak bulk interaction is added to the junction Hamiltonian, using the FRG-approach as we do in Sec. mu we 
still use Eq. (IC15I) for the charge- and the spin-conductance tensors by just replacing the ’’bare” S'-matrix elements 
with the running ones, Sjj'(D). For instance, this led us to write Eq. (13^1) of the main text. The result one obtains 
does actually correspond to the conductance tensor of the ’’connected” junction, that is, of the junction of weakly 
interacting wires connected to bulk Fermi liquid reservoirs. As pointed out in2 & 2LZ£, the tensor G c / S y t (j t ji) defined in 
Eq. (IC13I) are related to their analogs for the junction disconnected from the leads, G c (s)\(j.j') ■ by means of pertinent 
generalizations of Eq. (2.7) of Refs. [20ll2l| . that is. 


[G 1 ]c(s);(j,j') - [G 1 ] 


c(s);(Aj') + ^in,c(s)^d' 


(C16) 


with the interface charge (spin) conductance G; n c ( s ) = J 9c ^2i ~ ■ From Eq. (IC16I) one sees that the matrix elements 
G c ( s );(j,j') explicitly depend on the bulk interaction parameters via the Luttinger parameters g c ^ (see Fig. [9] for a 
sketch of the junction connected to reservoirs with generic values of the Luttinger parameter g R ). Equation (IC16I) 
suggests how to recover the conductance tensor for the junction disconnected from the leads within the FRG-approach: 
it is enough to just compute the G c i s \.^j ji)S and, therefore, to invert Eq. (IC16I) to derive the G^Wj^qs. 

We now discuss how Eq. (1C 161) is modified in the dual models we used to apply the FRG-approach to junctions at 
a strong bulk interaction in the wire. 


2. Junction of strongly interacting quantum wires at g c = g s = 3 

As Eq. (IG4I) applies independently of the specific form of the current operators, in order to generalize Eq. 1C16|) to 
the strongly interacting case we discuss in Sec. IIIIB 1 one has just to compute the 'D{j lCr y,(j>x' ,t')-terms 
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g L 


FIG. 9: Sketch of a connected junction of K wires with parameters g, u connected to reservoirs with parameters gL, ul- 

in Eq. (IC3I) using the current operators in the dual theory given in Eq. (l49l) . As a result, one now obtains 

with 

J dE f(E)[f(-E-u)- f(-E + u)\ 

x 1 2 [is.^.,_r |2 + |S i+1 , /+1 | 2 3 4 ] 

- e i ^-^[S j . w +S j+1 , r . 1 ] - e^(^[IVr^+rl 2 + |^ +1 , / _ 1 | 2 ]| . (C17) 

Following the same steps as in Appendix 1C 11 from Eq. (IC17I) one therefore readily derives the dc conductance tensor, 
which is given by 


G 




= S„.„,G 


0,3 


(CIS) 


with 


Gj,j' — ^ j 2 < ^b' + l^f+ij'+il 2 + j'-il 2 

- + $j- i.i'+i + I'Sj-iJ'+il 2 + l^+ij'-i| 2 ]| • (C19) 

From Eq. (IC19I) one therefore computes the conductance tensor at g c = g s = 3. When either g Cl or g s (or both) 
are different from 3, using Eq. (IC19I) one recovers the conductance tensor in the case in which g c and g s both 
asymptotically tend to 3 in the leads. The charge- and the spin-conductance tensors for the disconnected junction 
are obtained using Eq. (IC16I) with G; n . c ( s ) = v- 


1 L. D. Landau, Sov. Phys. JETP 3, 920 (1957). 

2 L. D. Landau, Sov. Phys. JETP 5, 101 (1957). 

3 L. D. Landau, Sov. Phys. JETP 8, 70 (1958). 

4 A. A. Abrikosov, L. P. Gorkov, and I. E. Dzyaloshinski, Methods of Quantum field theory in statistical physics (Dover, New 
York, 1963). 

M. P. A. Fisher and L. A. Glazman, in Mesoscopic Electron Transport , edited by G. S. L. Kouwenhoven and L. Sohn (NATO 
ASI Series E, Kluwer Ac. Publ., Dordrecht, 1997). 


5 














29 


6 P. Nozieres, Theory of interacting Fermi systems (Addison-Wesley, Reading (MA), 1964). 

7 J. Solyom, Adv. Phys. 28, 201 (1979). 

8 R. Shankar, Rev. Mod. Phys. 66, 129 (1994). 

9 J. Polchinski, in Proceedings of 1992 Theoretical Advanced Studies Institute in Elementary Particle Physics, edited by 
J. Harvey and J. Polchinski (World Scientific, Singapore, 1993). 

10 A. Furusaki and N. Nagaosa, Phys. Rev. B 47, 4631 (1993). 

11 C. L. Kane and M. P. A. Fisher, Phys. Rev. Lett. 68, 1220 (1992). 

12 C. L. Kane and M. P. A. Fisher, Phys. Rev. B 46, 15233 (1992). 

13 K. A. Matveev, D. Yue, and L. I. Glazman, Phys. Rev. Lett. 71, 3351 (1993). 

14 D. Yue, L. I. Glazman, and K. A. Matveev, Phys. Rev. B 49, 1966 (1994). 

15 Y. V. Nazarov and L. I. Glazman, Phys. Rev. Lett. 91, 126804 (2003). 

16 C. Nayak, M. P. A. Fisher, A. W. W. Ludwig, and H. H. Lin, Phys. Rev. B 59, 15694 (1999). 

17 S. Chen, B. Trauzettel, and R. Egger, Phys. Rev. Lett. 89, 226404 (2002). 

18 S. Lai, S. Rao, and D. Sen, Phys. Rev. B 66, 165327 (2002). 

19 K.-V. Pham, F. Piechon, K.-l. Imura, and P. Lederer, Phys. Rev. B 68, 205110 (2003). 

20 C. Chamon, M. Oshikawa, and I. Affleck, Phys. Rev. Lett. 91, 206403 (2003). 

21 M. Oshikawa, C. Chamon, and I. Affleck, J. Stat. Mech. P02008 (2006). 

22 S. Rao and D. Sen, Phys. Rev. B 70, 195115 (2004). 

23 C.-Y. Hou and C. Chamon, Phys. Rev. B 77, 155422 (2008). 

24 P. Jordan and E. Wigner, Z. Phys. 47, 631 (1928). 

25 S. Eggert and I. Affleck, Phys. Rev. B 46, 10866 (1992). 

26 I. Affleck, J. Phys. A 31, 2761 (1998). 

2 ‘ N. Crampe and A. Trombettoni, Nucl. Phys. B 871, 526 (2013). 

28 A. M. Tsvelik, Phys. Rev. Lett 110, 147202 (2013). 

29 S. A. Reyes and A. M. Tsvelik, Phys. Rev. Lett. 95, 186404 (2005). 

30 D. Giuliano, D. Rossini, P. Sodano, and A. Trombettoni, Phys. Rev. B 87, 035104 (2013). 

31 L. I. Glazman and A. I. Larkin, Phys. Rev. Lett. 79, 3736 (1997). 

32 D. Giuliano and P. Sodano, Nucl. Phys. B 837, 153 (2010). 

33 A. Cirillo, M. Mancini, D. Giuliano, and P. Sodano, Nucl. Phys. B 852, 235 (2011). 

34 B. Beri and N. R. Cooper, Phys. Rev. Lett. 109, 156803 (2012). 

35 M. R. Galpin, A. K. Mitchell, J. Temaismithi, D. E. Logan, B. Beri, and N. R. Cooper, Phys. Rev. B 89, 045143 (2014). 

36 B. Beri, Phys. Rev. Lett. 110, 216803 (2013). 

3 ‘ E. Eriksson, C. Mora, A. Zazunov, and R. Egger, Phys. Rev. Lett. 113, 076404 (2014). 

38 E. Eriksson, A. Nava, C. Mora, and R. Egger, Phys. Rev. B 90, 245417 (2014). 

39 S. Tomonaga, Prog. Theor. Phys. 13, 467 (1955). 

40 J. M. Luttinger, J. Math. Phys. 4, 1154 (1963). 

41 F. D. M. Haldane, J. Phys. C: Solid State Phys. 14, 2585 (1981). 

42 F. D. M. Haldane, Phys. Rev. Lett. 47, 1840 (1981). 

43 S. Das, S. Rao, and D. Sen, Phys. Rev. B 74, 045322 (2006). 

44 D. N. Aristov, Phys. Rev. B 83, 115446 (2011). 

45 P. Fendley, A. W. W. Ludwig, and H. Saleur, Phys. Rev. Lett. 75, 2196 (1995). 

46 P. Fendley, A. W. W. Ludwig, and H. Saleur, Phys. Rev. Lett. 74, 3005 (1995). 

47 P. Fendley, A. W. W. Ludwig, and H. Saleur, Phys. Rev. B 52, 8934 (1995). 

48 A. Rahmani, C.-Y. Hou, A. Feiguin, C. Chamon, and I. Affleck, Phys. Rev. Lett. 105, 226803 (2010). 

49 A. Rahmani, C.-Y. Hou, A. Feiguin, M. Oshikawa, C. Chamon, and I. Affleck, Phys. Rev. B 85, 045120 (2012). 

50 N. Sedlmayr, D. Morath, J. Sirker, S. Eggert, and I. Affleck, Phys. Rev. B 89, 045133 (2014). 

51 A. E. Feiguin and S. R. White, Phys. Rev. B 72, 220401(R) (2005). 

52 C. Karrasch, J. H. Bardarson, and J. E. Moore, Phys. Rev. Lett. 108, 227206 (2012). 

53 V. Meden, S. Andergassen, W. Metzner, U. Schollwock, and K. Schonhammer, Europhys. Lett. 64, 769 (2003). 

54 S. Andergassen, T. Enss, V. Meden, , W. Metzner, U. Schollwock, and K. Schonhammer, Phys. Rev. B 70, 075102 (2004). 

55 V. Meden, T. Enss, S. Andergassen, W. Metzner, and K. Schonhammer, Phys. Rev. B 71, 041302(R) (2005). 

56 T. Enss, V. Meden, S. Andergassen, X. Barnabe-Theriault, W. Metzner, and K. Schonhammer, Phys. Rev. B 71, 155401 
(2005). 

57 S. Andergassen, T. Enss, V. Meden, W. Metzner, U. Schollwock, and K. Schonhammer, Phys. Rev. B 73, 045125 (2006). 

58 X. Barnabe-Theriault, A. Sedeki, V. Meden, and K. Schonhammer, Phys. Rev. Lett. 94, 136405 (2005). 

59 V. Meden, S. Andergassen, T. Enss, H. Scholler, and K. Schonhammer, New Journal of Physics 10, 045012 (2006). 

60 X. Barnabe-Theriault, A. Sedeki, V. Meden, and K. Schonhammer, Phys. Rev. B 71, 205327 (2005). 

61 K. Janzen, V. Meden, and K. Schonhammer, Phys. Rev. B 74, 085301 (2006). 

62 S. G. Jakobs, V. Meden, H. Scholler, and T. Enss, Phys. Rev. B 75, 035126 (2007). 

63 M. Titov, M. Muller, and W. Belzig, Phys. Rev. Lett. 97, 237006 (2006). 

64 H. J. Schulz, Phys. Rev. B 22, 5274 (1980). 

65 D. N. Aristov and P. Woffle, Phys. Rev. B 86, 035137 (2012). 

66 D. N. Aristov and P. Woffle, Phys. Rev. B 88, 075131 (2013). 

67 D. N. Aristov and P. Woffle, Phys. Rev. B 84, 155426 (2011). 




30 


68 J. von Delft and H. Schoeller, Annalen Phys. 7, 225 (1998). 

69 C.-Y. Hou, A. Rahmani, A. E. Feiguin, and C. Chamon, Phys. Rev. B 86, 075451 (2012). 

70 S. Das, S. Rao, and A. Saha, Phys. Rev. B 77, 155418 (2008). 

71 I. Affleck and D. Giuliano, J. Stat. Mech. P06011 (2012). 

72 L. Fidkowski, J. Alicea, N. H. Lindner, R. M. Lutchyn, and M. P. A. Fisher, Phys. Rev. B 85, 245121 (2012). 

73 D. Giuliano and I. Affleck, J. Stat. Mech. P02034 (2013). 

74 D. Giuliano and I. Affleck, Phys. Rev. B 90, 045133 (2014). 

75 I. Affleck and D. Giuliano, J. Stat. Phys. , Special Issue in Memory of K. G. Wilson 157, 666 (2014). 

76 H. J. Shultz, G. Cuniberti, and P. Pieri, in Field Theory for Low-Dimensional Condensed Matter Systems , edited by V. T. 
G.Morandi, P. Sodano and A. Tagliacozzo (Springer-Verlag, 1997). 

77 A. O. Gogolin, A. A. Nersesyan, and A. M. Tsvelik, Bosonization and Strongly Correlated Systems (Cambridge University 
Press, 2004). 

78 T. Sugiyama, Prog. Theor. Phys. 64, 406 (1980). 

79 D. L. Maslov and M. Stone, Phys. Rev. B 52, R5539 (1995). 




